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ABSTRACT 


This  report  discusses  the  results  of  VHF  radio  direction  finding  (DF)  experinaents  carried 
out  with  an  experimental  linear  eight-antenna  DF  system.  The  main  purpose  of  these 
experiments  was  to  investigate  the  effect  of  the  environment,  in  particular  multipath  prop¬ 
agation,  on  DF  performance.  A  theoretical  analysis  of  all  known  error  mechanisms  has 
also  been  carried  out  in  order  to  better  understand  the  measured  results.  This  includes 
temporal  noise,  equipment  errors,  mutual  coupling,  and  multipath.  The  analysis  of  mul¬ 
tipath  also  includes  the  development  of  new  models  to  predict  statistical  effects.  The 
results  of  this  investigation  indicate  that  multipath  effects  are  complex  and  may  create 
difficulties  in  modem  systems  if  improved  accuracy  is  an  objective.  Further  measurements 
will  be  required  to  verify  these  preliminary  findings. 

RESUME 

Ce  rapport  traite  des  resultats  experimentaux  issus  de  tests  realises  a  I’aide  d’un  systeme 
experimental  de  radiogoniometrie  a  huit  antennes  VHF  repaxties  lineairement.  L’objectif 
principal  poursuivi  lors  de  ces  tests  etait  d’etudier  les  effets  environnementaux  sur  les 
performances  des  mesures  radiogoniometriques  et,  plus  paxticulierement,  I’effet  des  prop¬ 
agations  multi-voies.  Afin  d’obtenir  une  meilleure  comprehension  des  resultats  mesures, 
une  analyse  theorique  de  tons  les  types  connus  de  mechanismes  d’erreur  a  ete  menee, 
incluant  les  erreurs  dues  aux  bruits  variant  dans  le  temps,  aux  problemes  d’equipement, 
a  I’effet  de  couplage  mutuel  et,  enfin,  aux  propagations  multi-voies.  L’analyse  des  propa¬ 
gations  multi-voies  comprend  egalement  le  developpement  de  nouveaux  modeles  visant  a 
predire  les  effets  statistiques.  Les  resultats  de  ces  recherches  indiquent  que  les  effets  des 
propagations  multi-voies  s’averent  d’une  grande  complexite  et  causeront  probablement 
certains  problemes  aux  systemes  modernes  radiogoniometriques  cherchant  a  ameliorer  la 
precision  des  mesmes.  Ces  resultats  preliminaires  devront  toutefois  etre  verifies  pax  des 
tests  et  mesures  supplementaires. 
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EXECUTIVE  SUMMARY 


Improvement  of  communications  radio  direction  finding  (DF)  accuracy  is  a  high  priority 
for  the  Canadian  and  Allied  Forces.  To  this  end,  research  in  advanced  DF  techniques  has 
been  carried  out  worldwide  over  the  last  two  decades  with  the  view  of  taking  advantage 
of  advances  in  DF  algorithms  as  well  as  the  capabilities  of  modern  processing  technology. 
At  the  Defence  Research  Establishment  Ottawa  (DREO),  an  experimental  eight-channel 
VHF  DF  system,  called  the  Osprey  System,  has  been  set  up  to  investigate  the  distortion 
of  the  received  signal  wavefront  caused  by  multipath  propagation.  Field  trials  with  this 
system  took  place  in  the  Fall  of  1992  and  the  Spring  of  1993.  These  trials  were  carried 
out  in  support  of  the  Canadian  Forces  for  the  improvement  of  land  tactical  VHF /UHF 
DF  accuracy,  and  in  support  of  TTCP  QTP-18  to  investigate  environmental  effects  on  a 
VHF  DF  system. 

This  report  discusses  and  analyzes  the  results  of  these  field  trials.  The  theoretical 
effects  of  error  mechanisms  which  degraded  the  DF  measurements  are  analyzed  in  detail, 
including  the  modelling  of  the  multipath  environment.  This  analysis  is  then  used  to 
break  down  the  measured  errors  in  terms  of  these  mechanisms.  The  effects  due  to  each 
mechanism  are  then  studied  in  more  detail. 

Quantitatively,  multipath  propagation  was  found  to  be  the  greatest  source  of 
error.  It  was  followed  closely  by  mutual  coupling.  Based  on  the  modelling  studies  of  both 
phenomena,  important  qualitative  differences  were  found.  Mutual  coupling  errors  are  a 
function  of  transmitter  azimuth  only  (when  a  plane  wave  signal  is  considered)  and  are 
predictable.  The  error  pattern  changes  relatively  slowly;  cycling  once  from  the  maximum 
positive  error  to  the  maximum  negative  error  over  a  range  of  180  degrees  in  azimuth. 
Multipath  errors  are  a  function  of  both  transmitter  azimuth  and  range.  The  error  pattern 
becomes  more  unpredictable  (noise-like)  as  the  range  of  the  multipath  source  from  the  DF 
system  increases  (at  a  range  of  0  meters  multipath  and  mutual  coupling  can  be  regarded  as 
the  same  phenomena).  These  qualitative  differences  make  it  relatively  easy  to  distinguish 
mutual  coupling  effects  from  multipath  effects. 

In  examining  the  actual  multipath  induced  errors,  two  important  observations 
were  made.  The  first  was  that  the  multipath  errors  were  very  noise-like  indicating  sig¬ 
nificant  sources  of  multipath  at  ranges  of  several  hundred  meters  or  more.  The  second 
was  that  there  appears  to  be  a  transmitter  range  dependency  even  for  transmitter  ranges 
over  5  km.  Based  on  the  theoretical  models,  this  implies  that  sources  of  multipath  which 
are  reflective  in  nature  (i.e.  redirect  the  incident  signal  in  a  new  direction  with  limited 
scattering)  are  significant.  Both  observations  lead  to  complications  for  site  calibration 
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if  improved  DF  accuracy  is  the  objective.  They  indicate  that  calibration  of  the  array 
steering  vector  (or  array  manifold)  at  a  particular  frequency  needs  to  be  performed  as  a 
function  of  both  range  and  azimuth  ,  not  azimuth  as  is  usually  done.  Further  experimen¬ 
tation  will  be  required  to  verify  these  observations  before  any  definite  conclusions  can  be 
drawn. 

In  general  the  results  described  in  this  report  indicate  that  the  multipath  envi¬ 
ronment  is  very  complex.  Modelling  can  be  used,  with  some  success,  to  better  understand 
the  statistical  effects  of  multipath.  For  future  work,  the  theoretical  models  need  to  be 
refined  to  better  approximate  real  life  effects.  Further  field  testing  will  also  be  required 
to  verify  these  models  and  determine  appropriate  model  parameters  (e.g.  density  of  mul¬ 
tipath  sources,  relative  importance  of  reflective  versus  scattering  types,  etc.).  This  could 
ultimately  lead  to  the  development  of  new  approaches  to  VHF  DF  including  (but  not 
restricted  to)  new  algorithms  and  new  calibration  techniques. 
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1.0  INTRODUCTION 


Improvement  of  communications  radio  direction  finding  (DF)  accuracy  is  a  high  priority 
for  the  Canadian  and  Allied  Forces.  To  this  end,  research  in  advanced  DF  techniques  has 
been  carried  out  worldwide  over  the  last  two  decades  with  the  view  of  taking  advantage 
of  advances  in  DF  algorithms  as  well  as  the  capabilities  of  modern  processing  technology. 
Central  to  this  approach  is  the  N-channel  digital  beamformer,  one  possible  version  of 
which  is  shown  in  block  diagram  form  in  Figure  1.  The  main  advantage  of  this  approach 
is  that  the  phase  and  amplitude  measurements  from  each  antenna  are  available  for  anal¬ 
ysis.  This  maximizes  the  information  available  about  the  incoming  radio  signal  which 
allows  multiple  signal  DF  (superresolution)  to  be  performed,  or  distortion  of  the  received 
wavefront  (compared  to  theoretical  expectations)  to  be  measured. 


Figure  1:  Experimental  N  channel  DF  system  block  diagram. 


At  DREO  an  eight-channel  hardware  realization  of  Figure  1,  called  the  Osprey  System, 
has  been  set  up  to  investigate  the  distortion  of  the  received  signal  wavefront  caused 
by  multipath  propagation.  This  report  describes  and  analyzes  the  results  of  VHF  DF 
field  trials  which  took  place  during  the  Fall  of  1992  and  the  Spring  of  1993.  The  trials 
were  carried  out  in  support  of  the  Canadian  Forces  for  the  improvement  of  land  tactical 
VHF/UHF  DF  accuracy,  and  in  support  of  TTCP  QTP-18  to  investigate  environmental 
effects  on  a  VHF  DF  system. 

This  report  is  arranged  in  six  main  sections.  The  first  section  is  this  introduction. 
Section  2  describes  the  experimental  setup  including  a  description  of  the  site  surrounding 
the  Osprey  System,  the  transmit  and  receive  hardware,  and  the  measurement  test  plan. 
Section  3  describes  the  processing  of  the  data  from  conversion  to  IQ  format  to  the  cal¬ 
culation  of  the  transmitter  signal  bearing.  Section  4  provides  a  theoretical  evaluation  of 
the  sources  of  error  including  the  development  of  multipath  models  for  the  assessment 
of  the  effects  of  multipath  on  DF.  Section  5  describes  the  results  in  terms  of  the  error 
sources  with  the  objective  of  quantifying  and  qualifying  the  effects  of  multipath.  Finally 
the  conclusions  are  presented  in  Section  6. 

To  avoid  ambiguities,  the  following  definitions  for  antenna  systems  have  been  used  in 
this  report: 

Antenna  A  structure  which  collects  electromagnetic  energy  and  converts  it 
into  electrical  energy.  It  has  a  single  output  which  is  connected  to  a  single 
receiver  channel. 

Antenna  Element  A  part  of  the  antenna  structure,  e.g.  a  vertical  dipole 
antenna  has  an  upper  element  and  a  lower  element. 

Antenna  Array  A  collection  of  N  antennas  which  are  connected  to  N  re¬ 
ceiver  channels. 

Note  that  these  definitions  do  not  necessarily  reflect  the  most  common  usuage  of  these 
terms.  They  have  been  defined  in  a  way  that  is  most  appropriate  for  this  report. 

2.0  EXPERIMENTAL  SETUP 

2.1  Site  Layout 

The  layout  of  the  site  was  as  shown  in  Figure  2.  This  was  neither  flat  nor  cleared  of 
obstacles,  but  rather  a  site  located  on  a  flat  field  at  DREO  surrounded  by  trees,  low  wire 
fences,  roads,  a  ditch,  etc.  In  short,  a  variety  of  obstacles  that  according  to  [1]  can  cause 
errors  in  bearing  measurements  but  are  typical  of  the  Ottawa  area. 
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Signal 


Figure  2:  Physical  layout  of  field  trial  DF  site. 

The  Osprey  antenna  system  consisted  of  a  linear  array  of  2.9  meter  dipole  whip  an¬ 
tennas  spaced  at  intervals  of  2  meters.  The  antennas  were  mounted  on  a  15  meter  long 
fiberglass  truss.  The  truss  was  in  turn  supported  2  meters  off  the  ground  by  a  large  pivot 
at  the  center,  enabling  the  truss  to  be  rotated  manually  through  360  degrees.  Adjustable 
supports  at  either  end  of  the  truss  were  used  to  prevent  sagging  over  the  length  of  the 
truss.  For  the  experiments  discussed  in  this  report,  the  truss  was  left  in  the  position 
shown  in  Figure  2  and  the  0  degree  azimuth  bearing  reference  was  also  as  shown. 

For  alignment  purposes,  the  0  degree  bearing  reference  was  chosen  to  correspond  to 
the  bearing  of  a  water  tower  located  at  a  range  of  990  meters  from  the  antenna  system 
(not  shown  in  Figure  2)  which  provided  a  useful  map  and  visual  marker.  The  alignment 
accuracy  was  assessed  to  be  ±1  degree. 

2.2  Acquisition  System 

The  data  acquisition  system  shown  in  Figure  1  was  designed  for  operation  in  the  30-100 
MHz  band.  Each  receiver  consists  of  two  parts:  2m  RF  front  end  mounted  at  the  base  of 
the  antenna  to  downconvert  the  signal  to  10.7  MHz,  and  a  rack  mounted  IF  unit  located 
in  the  acquisition  trailer  to  downconvert  the  10.7  MHz  signal  to  baseband.  A  common 
reference  signal  was  used  to  control  the  input  RF  frequency  of  the  front  end  units.  Given 
its  design,  there  were  actually  two  input  frequency  bands,  namely,  frej  —  10.7  MHz  and 
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frej  + 10.7  MHz  where  fref  was  the  frequency  of  the  reference  signal.  For  the  experiments, 
the  reference  frequencies  /,*/  =  70  MHz  and  fref  =  51.8  MHz  were  used  in  the  Fall  and 
Spring  respectively.  To  eliminate  one  of  the  input  frequency  bands,  55-67  MHz  band  pass 
filters  were  used  at  the  input  of  the  RF  units.  These  filters  also  eliminate  strong  signals 
in  the  FM  band  which  would  have  otherwise  resulted  in  the  generation  of  a  high  level  of 
intermodulation  noise  in  the  receivers. 

A  digitally  controlled  attenuator  in  each  receiver  was  used  to  adjust  the  IF  gain.  The 
attenuator  settings  were  controlled  by  the  A/D  which  was  in  turn  controlled  by  a  486  PC. 
During  the  experiments  the  attenuators  were  left  in  the  0  dB  setting,  since  changing  the 
attenuator  settings  was  found  to  affect  the  relative  phase  offsets  between  channels. 

Receiver  tuning  was  controlled  through  a  single  common  synthesizer  which  is  currently 
manually  controlled  but  can  be  modified  for  computer  control.  The  baseband  output  of  the 
receivers  is  low  pass  filtered  with  an  effective  noise  bandwidth  of  7.1  kHz.  This  output 
is  fed  to  the  A/D  converter  which  synchronously  samples  from  one  to  eight  channels. 
The  digital  output  data  is  stored  on  the  PC  where  it  can  be  further  processed.  The 
channels  sampled,  sampling  rate,  and  data  blocksize  are  all  computer  controlled.  For  the 
experiments  all  eight  channels  were  recorded  at  a  sampling  rate  of  1.0  kHz.  The  blocksize 
was  also  limited  to  2048  samples  for  each  measurement  to  minimize  the  amount  of  data 
that  required  storage. 

For  the  field  trials,  the  transmitter  equipment  was  placed  in  a  Dodge  minivan  and  a  2 
meter  monopole  transmitting  antenna  attached  to  the  roof  of  the  vehicle.  The  transmitter 
equipment  consisted  of  a  synthesizer  and  power  amplifier  capable  of  generating  a  CW 
signal  of  up  to  3  watts  of  power. 

2.3  The  Test  Plan 

The  field  trials  consisted  of  four  days  of  data  collection  which  were  the  and  6‘**  of 
November  1992  and  the  24‘^  and  25**^  of  March  1993.  The  routes  followed  on  these  days 
are  shown  in  Figures  3-6.  For  each  trial  the  transmitter  was  positioned  at  the  start  of 
the  route  imtil  a  DF  measurement  (or  measurements)  was  taken,  then  moved  50  meters 
along  the  route  to  the  next  position  for  the  next  measurement.  This  was  repeated  until 
the  whole  route  was  covered.  Moving  the  transmitter  and  taking  a  DF  measurement 
generally  took  about  30  seconds. 

For  logistical  reasons  the  chosen  routes  followed  local  roads,  and  the  transmitter  was 
aiway  positioned  on  the  right  side  of  the  road  relative  to  the  direction  of  travel.  In  a 
number  of  cases  the  routes,  or  parts  of  the  routes,  were  repeated  on  the  same  or  following 
day  to  investigate  changes  in  the  measurements  over  time. 
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The  start  position  for  each  route  was  chosen  to  be  easy  to  locate  on  a  1:25000  to¬ 
pographical  map  of  the  Ottawa  area  or  a  1:5000  overhead  photograph  of  the  DREO 
site.  Distance  from  the  start  position  was  measured  using  a  fifth  wheel  (a  bicycle  wheel 
mounted  on  the  back  bumper  of  the  vehicle  which  measures  distance  travelled)  whose 
smallest  measurement  increment  was  10  meters. 

On  the  5‘**  of  November  and  the  24***  of  March,  one  DF  measurement  was  taken  for 
each  transmitter  position.  On  the  6***  of  November  two  consecutive  measurements  were 
taken  at  a  time.  The  objective  was  to  determine  the  repeatability  of  the  measurements 
for  each  transmitter  position.  Unfortunately,  due  to  synchronization  problems  the  trans¬ 
mitter  was  sometimes  moving  during  the  second  measurement.  When  this  happened  the 
second  measurement  was  discarded.  On  the  25***  of  March  three  consecutive  DF  mea¬ 
surements  were  taken  at  a  time.  Again,  due  to  synchronization  problems,  some  of  these 
measurements  had  to  be  discarded. 

The  transmitter  frequencies  used  were  59.3  MHz  during  the  November  trials,  and  62.5 
MHz  during  the  March  trials.  Slight  adjustments  to  these  frequencies  were  made  to  ensure 
the  final  baseband  frequency  of  the  received  signal  of  the  Osprey  System  was  between  150 
Hz  and  500  Hz. 

During  the  November  trials  the  temperatures  were  relatively  constant  and  remained 
at  about  2°  C,  the  sky  was  overcast,  and  the  ground  was  wet.  During  the  March  trials 
the  temperatures  ranged  from  0°  C  in  the  morning  to  +7“  C  in  the  afternoon,  the  sky 
was  clear,  and  the  ground  was  frozen  and  snow  covered  to  a  depth  of  50  cm. 
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Figure  3:  Transmitter  locations  (marked  by  +)  on  November  5,  1992  relative  to  Os¬ 
prey  system  at  (0,0).  Arrows  show  vehicle  direction.  Start  position  was  (xE-w,yN-s)  = 
(-5825,1909) 
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Figure  4:  Transmitter  locations  (marked  by  +)  on  November  6,  1992  relative  to  Osprey 
system.  Arrows  show  vehicle  direction.  Start  position  was  (xE-WjyN-s)  =  (—5825, 1909) 
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Figure  5:  Transmitter  locations  (marked  by  +)  on  March  24,  1993  relative  to  Osprey 
system.  Arrows  show  vehicle  direction.  Start  position  was  (xe-WiYn-s)  =  (45,-8) 
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Figure  6:  Treinsmitter  locations  (marked  by  +)  on  March  25,  1993  relative  to  Osprey 
system.  Arrows  show  vehicle  direction.  Start  position  was  (xE-w,yN-s)  =  (118,110) 


3.0  SIGNAL  PROCESSING  CONSIDERATIONS 


3.1  IQ  Conversion 

An  example  of  the  positive  frequency  spectrum  of  the  raw  data  is  shown  in  Figure  7. 
Here  a  2048  point  FFT  (with  a  Blackman  window)  was  performed  on  the  raw  data  (2048 
samples)  collected  from  antenna  0.  The  received  (down  converted)  CW  signal  is  clearly 
evident  at  322  Hz  along  with  some  interfering  signals  (likely  receiver  generated)  and  the 
received  noise. 


For  the  DF  algorithm  discussed  in  this  report  (Section  3.3),  it  is  more  convenient  to 
have  the  data  in  complex  (I  and  Q)  baseband  format.  Although  many  different  approaches 
for  accomplishing  this  are  possible,  the  approach  used  here  is  to  take  advantage  of  the 
fact  that  the  signal  used  was  CW.  This  allows  the  use  of  very  narrow  band  filters  centered 
on  the  signal.  Noise  and  interference  outside  the  filter  passband  are  therefore  rejected, 
significantly  improving  the  signal-to-noise  ratio  (SNR). 
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The  actual  approach  is  based  on  the  discrete  Fourier  transform.  The  definition  of  the 
discrete  Fourier  transform  of  a  data  sequence  x,(fc)  for  k  =  0,l,...,K  —  1  is  given  by, 

S,(U)  =  E  (1) 

Jh=0 

where  t  is  the  antenna  number  (0  to  —  1),  /„  is  the  discrete  frequency  such  that 

/«  =  ^  =  (2) 

and  f,  is  the  sampling  frequency.  The  transformed  data  represents  the  I  and  Q  values 
of  the  original  data  measured  at  all  possible  values  of  /«.  For  a  CW  signal,  it  is  only 
necessary  to  determine  S{fo)  where  fo  is  the  frequency  of  the  CW  signal. 

One  difficulty  is  that  the  frequency  spectrum  can  only  be  calculated  at  discrete  points 
which  are  multiples  of  fa/K.  To  improve  this  situation,  zeros  padding  can  be  used  (i.e. 
appending  zeros  to  the  end  of  the  data)  to  increase  the  number  of  sample  points  from  K  to 
K'.  In  the  limit  where  K'  —*  oo,  the  frequency  becomes  a  continuous  function  and  Si{f) 
can  be  c<dculated  for  any  frequency  /  such  that  0  <  f  <  fa-  Note  that  since  Xi{k)  =  0  for 
k>  K,  the  upper  limit  of  the  summation  in  equation  (1)  will  still  be  —  1,  not  K'  —  1. 

In  general,  the  exact  frequency  of  the  CW  signal  will  be  unknown  so  that  the  first 
step  is  to  perform  the  Fourier  transform  without  zero  padding  {K'  =  K).  Using  the 
Fast  Fourier  Transform  (FFT)  to  perform  the  computations,  the  approximate  frequency 
fo  is  taken  as  the  frequency  f  =  f^  which  maximizes  the  value  of  |5,(/„)|.  Using  this 
approximate  value  as  a  starting  point,  equation  (1)  is  then  used  to  fine  time  the  results. 

The  above  approach  assumes  that  the  signal-to-interference-plus-noise  ratio  (SNIR)  is 
high  (as  in  Figure  7)  in  all  of  the  receive  channels,  which  was  true  for  the  experiments 
described  in  this  report.  The  choice  of  which  channel  i  to  use  is  arbitrary,  although 
for  consistency,  genereilly  i  =  0  was  used.  At  lower  signal  levels  a  more  accurate  (but 
computationally  more  expensive)  approach  would  be  to  maximize  the  sum  of  the  power 
spectra,  that  is,  maximize  the  function 

E  (3) 

»=o 

where  N  is  the  number  of  sensor  channels. 

Once  the  data  has  been  converted  to  the  complex  baseband  format,  a  single  sample 
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or  snapshot  of  all  N  channels  is  represented  by 


(4) 


Xn-1 


where  Xi  =  5,(/o)  for  i  =  0, —  1  (see  also  equation  (1)).  Note  that  in  the  IQ 
processing  scheme  described  above  K  real  data  points  per  channel  are  collected  to  generate 
a  single  snapshot  in  IQ  form. 

For  consecutive  measurements  measured  at  times  starting  at  to,  ti,  t2,  ...,tT~u  a  data 
matrix  can  be  formed  as, 


X  =  [x(to),  x(ti),  x(t2), . .  • ,  x(tr_i)]  (5) 

where  x(tn)  represents  a  snapshot  measured  at  time  and  T  is  the  number  of  snapshots. 
3.2  Channel  Equalization 

Channel  equalization  is  normally  performed  in  direction  finding  to  ensure  that  the  channel 
gains  and  phases  are  matched.  This  is  often  done  by  disconnecting  the  receive  antennas 
and  replacing  the  antenna  feeds  by  test  signals  of  identical  frequencies,  amplitudes,  and 
phases.  Adjustable  gain  and  phase  devices  are  then  used  to  adjust  the  channels  until  the 
measured  receiver  outputs  are  identical  in  both  gain  and  phase.  This  procedure,  however, 
does  not  account  for  mismatches  introduced  by  the  antennas  themselves. 

An  alternate  approach,  which  has  been  adopted  in  this  report,  is  to  assume  the  channel 
gain  and  phase  mismatches  introduced  by  the  equipment  were  unknown,  and  then  estimate 
these  values  based  on  the  measured  data  collected  and  some  knowledge  about  the  signal 
environment.  Accurate  knowledge  about  the  signal  environment  is  critical,  otherwise 
errors  are  introduced  into  this  approach. 

Channel  equalization  can  be  performed  on  the  measured  data  X  using  the  expression 

Y  =  A-*X  (6) 

where  Y  is  the  corresponding  N  x  T  matrix  of  equalized  data,  and  A  is  an  N  x  N 
diagonal  matrix  with  complex  gain  elements  cq,  oj,  02, . . . ,  as-i-  These  complex  elements 
can  also  be  represented  in  terms  of  the  real  valued  gain  and  phase,  i.e.  for  channel  k 
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(0<k<N-l) 


ak  =  5*6®*  (7) 

where  5*  is  the  real-valued  gjun  and  6k  the  phase.  Since  5*  and  6k  are  both  relative  terms, 
this  implies  c*  will  be  ambiguous  unless  suitable  references  or  constraints  are  chosen.  To 
remedy  this  problem  the  phase  is  referenced  to  channel  0  so  that 

=  0  (8) 


and  the  gains  au-e  normalized  so  that 

N~1 

Tsk^N  (9) 

*=o 

Determination  of  the  complex  gain  elements  co,  Cj, . . . ,  otv-i  is  relatively  straight  for¬ 
ward  if  both  X  and  Y  are  known.  In  reality  Y  will  not  be  completely  known  and  therefore 
must  be  estimated  instead.  This  estimation  requires  a  complete  knowledge  of  the  antenna 
array  and  the  environment  in  which  it  is  working.  This  includes  the  effects  of  temporal 
noise,  mutual  coupling,  and  the  effects  of  multipath.  Temporal  noise  can  be  minimized 
by  performing  the  measurements  at  high  SNR’s  (>  30  dB).  Mutual  coupling  effects  can 
be  measured  or  modelled,  however,  they  are  ignored  in  this  analysis;  the  consequences  of 
which  are  discussed  in  Section  4.4.  Multipath  effects  are  not  suflBciently  well  understood 
that  they  can  be  measured  or  modelled,  hence  the  estimation  of  Y  is  problematic.  How¬ 
ever  in  this  analysis  special  precautions  have  been  taken  so  that  multipath  effects  can  be 
adequately  represented  by  spatial  white  noise. 

It  follows  from  these  assumptions  that  the  equalized  array  response  can  be  decomposed 
into  a  signal  and  a  noise  component  as 

Y  =  A-*X  =  SB  -I-  N  (10) 

where  S  is  the  N  xT  signal  matrix  which  is  described  in  more  detail  below,  B  is  a  diagonal 
TxT  matrix  of  complex  signal  amplitues,  also  described  below,  and  N  is  the  N  xT  noise 
matrix  which  represents  the  error  due  to  temporal  and  spatial  noise  effects. 

The  signal  matrix  S  can  be  defined  in  terms  of  its  columns  as 

S  =  (so,Si,S2,...,8r_i]  (11) 

where  So,Si,...,Sr_i  axe  N  x  1  signal  vectors  normalized  so  that  |sjt|  =  1  and  the  first 
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element  Sko  is  real- valued  for  0  <  ib  <  T  —  1.  Each  signal  vector  represents  the  expected 
response  of  the  antenna  array  to  the  signals  measured  during  the  corresponding  snapshot 
assuming  the  array  is  equalized  and  there  is  no  noise.  If  only  a  single  signal  is  being 
measured  during  each  snapshot,  and  the  antennas  are  omnidirectional,  then 


Sk  = 


Vn 


* 


gin^* 


■  i(N-t)frd 


tm4>k 


(12) 


where  d  represents  the  spacing  between  antennas  of  the  DF  array,  and  hd  sin  <l>k  represents 
the  path  length  difference  between  sensor  n  and  sensor  0.  The  assumptions  have  also  been 
made  that  the  signal  is  a  ground  wave  with  an  elevation  angle  of  0°,  and  that  the  shape  of 
the  wavefront  is  planar  (i.e.  the  transmitter  is  in  the  farfield  of  the  DF  array,  and  mutual 
coupling  and  multipath  effects  are  ignored).  In  the  special  case  of  a  single  signal,  is 
also  called  an  array  steering  vector. 

The  matrix  B  is  a  diagonal  matrix  with  diagonal  elements  6o,  •••,  hj-i-  The  element 

hk  represents  the  complex  amplitude  of  the  signal  vector  s*  in  the  snapshot. 

The  solution  of  equation  (10)  in  terms  of  the  complex  gains  A  involves  the  simultaneous 
solution  of  the  complex  signal  amplitudes  in  B,  since  neither  the  signal  amplitude  or 
phase  will  be  known  at  the  receiving  site.  The  method  by  which  this  equation  is  solved 
is  discussed  in  Appendix  A  and  is  called  the  stochastic  calibration  method  in  this  report. 
Since  this  is  a  statistical  approach  (least  squares  estimation)  the  greater  the  number  of 
measurements  the  better  the  estimate  of  the  channel  gains  and  phases.  This  also  includes 
varying  the  measurements  over  a  wide  range  of  azimuth  angles  if  spatial  noise  effects  are 
to  be  minimized. 


3.3  The  MUSIC  DF  Algorithm 

The  Multiple  Signal  Classification  (MUSIC)  DF  algorithm  has  been  extensively  reported 
in  the  open  literature  (for  example  see  references  [5]-[7])  and  is  a  popular  choice  among 
current  superresolution  algorithms  due  to  its  performance.  It  is  only  briefly  described 
here. 

Central  to  the  MUSIC  algorithm  is  the  N  x  N  covariance  matrix  defined  by 

R  =  E{yy"}  (13) 
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where  y  represents  a  single  snapshot  of  equalized  data  (i.e.  a  single  column  of  Y),  and 
the  superscript  ^  denotes  the  conjugate-transpose  operation.  Generally  this  quantity  is 
not  known  exactly  but  instead  estimated  from  the  data  using 

R  =  -^YY"  (14) 

M 

The  covariance  matrix  is  then  decomposed  in  terms  of  its  eigenvectors  and  eigenvalues  as 


N-l 

t=0 


where  A,-  represents  one  of  the  eigenvalues  (which  are  ordered  so  that  Aq  >  Ai  >  . . .  >  As), 
and  Vj  represents  the  corresponding  eigenvector  which  are  normalized  so  that  v^v,-  =  1. 

If  it  is  known  that  M  signals  are  being  received  at  the  tuned  frequency  then  the 
eigenvectors  can  be  divided  into  the  so-called  signal  subspace  eigenvectors  (vq,  Vjif_i) 
and  the  noise  subspace  eigenvectors  (va/,  v//_i).  The  MUSIC  spectrum  is  computed 

using  either  the  noise  subspace  eigenvectors  as 


_ 1 _ 


(16) 


or  the  signed  subspace  eigenvectors  as 


5(«  = 


1 


(17) 


Here,  S{^)  is  the  value  of  the  spectrum  at  the  given  azimuth  angle,  and  e(^)  is  the  steering 
vector  defined  earlier  in  Section  3.2  (e.g.  see  equation  (12)).  The  M  largest  peaks  in  the 
computed  spectrum  are  then  taken  to  correspond  to  the  direction  of  arrival  of  the  M 
signals. 


4.0  SOURCES  OF  ERROR 


4.1  Transmitter  Position 

The  determination  of  the  bearing  errors  in  the  DF  estimates  is  based  on  comparing  the 
measured  bearing  with  the  actual  bearing.  The  actual  bearing  can  be  determined  from 
the  relative  position  of  the  transmitter  with  respect  to  the  receiver.  Defining  the  position 
of  the  Osprey  elements  in  Cartesian  coordinates,  the  position  of  antenna  nfoi  0  <n  <  N 
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is  given  by 


(Xn,  yn)  =  {nd- i(A?  -  l)d,  0) 


(18) 


Using  the  center  of  the  antenna  array  (x,y)  =  (0,0)  as  a  reference,  the  bearing  of  a 
transmitter  at  the  location  {xtx,ytx)  is  easily  found  using 


<f>  =  arctan 


fe) 


(19) 


where  the  bearing  is  measured  in  a  clockwise  direction  from  the  y-axis  and  the  range  is 
found  using 

r  =  +  y<x  (20) 

,  Note  that  to  distinguish  position  parameters  of  the  form  “x”  and  “y”  from  data  parameters 
of  the  form  “i”  and  “y”,  only  the  data  parameters  are  italicized. 

In  real  life,  the  positions  of  transmit  and  receiving  antennas  will  not  be  known  exactly. 
This  introduces  errors  into  the  determination  of  the  actual  bearing.  The  amount  of  error 
is  dependent  on  the  manner  in  which  the  position  measurements  are  made. 

In  the  experiments  described  in  this  report,  the  transmitter  was  moved  to  successive 
positions  along  a  roadway.  The  equipment  used  to  measure  the  distance  between  these 
positions  had  a  minimum  resolution  of  10  m.  This  gives  a  position  accuracy  of  ±5  m  in 
the  direction  of  the  roadway.  Position  errors  due  to  other  sources  were  assessed  to  be 
small  enough  that  they  could  be  ignored.  Under  these  conditions  the  position  error  can 
be  considered  to  be  uniformly  distributed  over  the  ±5  m  interval  with  a  mean  of  0  m. 
The  corresponding  RMS  position  error  is  pt  =  10/\/T2  =  2.89  m.  Taking  the  worst  case 
(when  the  roadway  was  perpendicular  to  the  signal  direction),  the  RMS  bearing  error  can 
be  calculated  using 

=  arctan  (21) 

assuming  that  r  ^ 

To  avoid  confusion  in  the  rest  of  this  report,  the  term  “actual  bearing”  is  used  when  it 
is  assumed  the  transmitter  bearing  is  known  exactly  (i.e.  theoretical  derivations),  and  the 
term  “map  bearing”  is  used  when  only  an  estimate  is  known  (i.e.  experimental  results). 
Hence  p^  represents  the  RMS  map  bearing  estimation  error. 


4.2  Noise 

The  effect  of  noise  (temporal  or  spatial)  is  to  corrupt  the  phase  and  amplitude  measure¬ 
ments  made  at  each  antenna,  which  in  turn  degrades  the  accuracy  of  the  DF  estimation. 
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An  analysis  of  the  optimum  estimation  error  for  a  linear  array  against  a  single  signal  is 
given  in  [5].  Assuming  the  sensor  noise  is  additive  white  Gaussian  in  nature,  the  mean 
error  in  the  spatial  frequency  estimates  is  zero  and  the  variance  given  by, 


snrNiN^  - 1) 
where  the  spatial  frequency,  cj,,  is  defined  here  as 

27rd  . 

ijJt  =  — r-  sm  q> 

A 


(22) 


(23) 


and  snr  is  the  signal  to  noise  power  ratio.  In  terms  of  a  single  equalized  snapshot 

y  =  6s  +  n  (24) 

the  signal  to  noise  power  ratio  can  be  defined  as 

|6p 


snr  = 


n 


(25) 


where  6  is  the  complex  signal  amplitude,  s  is  the  signal  vector  defined  by  equation  (12), 
and  n  is  the  x  1  noise  vector. 

In  equation  (22),  the  variance  is  constant  with  respect  to  bearing.  To  determine  the 
corresponding  mean  error,  and  variance,  cr^,  of  the  bearing  estimates,  the  expressions 


and 


<f>e=  I  fu,{T)<f>e{T)dT 

<—00 


(26) 


(27) 


can  be  used  where  /u,(r)  is  the  probability  distribution  function  of  the  spatial  frequency 
errors,  and  <i>e{T)  is  the  bearing  error  given  by 


<1  (28) 
>1  (29) 


liquation  (29)  is  a  result  of  the  fact  that  the  bearings  for  a  linear  array  are  real- valued  and 
restricted  to  the  range  ±90”.  For  example,  given  a  transmitter  be2iring  of  89°  and  a  low 


(2^<“‘ 

for 

Vk*  +  ’■ 

for 
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value  of  snr,  the  maximum  positive  bearing  error  cannot  exceed  1°  while  the  maximum 
negative  bearing  error  could  be  as  much  as  —179°.  The  bearing  errors  in  this  case  are 
very  as3Tnmetric  with  respect  to  zero  error. 

Assuming  snr  >  1,  then  /^(t)  can  be  reasonably  approximated  as  having  the  same 
distribution  as  the  input  noise,  namely,  a  zero-mean  white  Gaussian  distribution.  There¬ 
fore, 

Ur)  =  (30) 

Since  this  leads  to  integrations  (equations  (26)  and  (27))  which  are  difficult  to  perform 
algebraically,  they  were  performed  numerically  and  the  results  shown  in  Figure  8.  In  this 
figure  the  absolute  value  of  the  mean  errors  (|<^e|)  and  the  standard  deviation  (a^)  of 
the  bearing  errors  measured  (theoretically)  using  the  Osprey  array  at  /  =  62.5  MHz  are 
plotted  as  a  function  of  signal  bearing  for  several  different  signal  to  noise  ratios.  Noise 
was  assumed  to  be  the  only  source  of  error. 

The  theoretical  RMS  bearing  error  is  given  by 

<l>eRMS  =  \l +  <^4,  (31) 

From  the  results  shown  in  Figure  8,  the  mean  error  is  significantly  smaller  than  the 
standard  deviation  for  the  same  bearing  and  SNR.  therefore  a  reasonable  approximation 
for  the  RMS  bearing  error  is  given  by 


4>eRMS  =  <^it>  (32) 

From  these  results  it  is  also  clear  that  the  best  accuracy  (lowest  RMS  error)  is  achieved 
for  the  broadside  direction  (0°  azimuth)  of  the  array  and  that  accuracy  significantly 
degrades  towards  the  endfire  directions  (±90°  azimuth).  The  “corner  effect”  that  occurs 
near  endfire  (e.g.  at  approx  75°  on  the  SNR  =  0  dB  curve)  is  due  to  the  asymmetry  effect 
described  earlier. 

For  conditions  where  the  signal  to  noise  ratio  is  greater  than  0  dB  (snr  >  1),  the 
accuracy  of  the  Osprey  array  can  be  reasonably  approximated  by  simpler  expressions. 
Away  from  the  endfire  directions,  a  reasonable  approximation  of  the  accuracy  is  given  by. 

In  the  endfire  position,  the  term  cos  <f)  goes  to  0  and  the  above  expression  for  accuracy  no 
longer  holds.  Instead,  it  has  been  found  empiracally  that  accuracy  can  be  approximated 
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Mean  Error  (degrees)  Standard  Deviation  (degrees) 


(b) 


(a) 

Figure  8:  Theoretical  absolute  mean  bearing  error  (a)  and  bearing  variance  (b)  of  the 
Osprey  System  as  a  function  of  signal  bearing  and  SNR. 
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17*  = 


(34) 


It  is  interesting  to  note  that  in  the  endfire  case,  the  accuracy  is  proportional  to  snr~'i 
compared  to  snr~i  for  the  broadside  case. 

Although  the  analysis  of  noise  effects  has  been  carried  out  using  white  Gaussian  noise 
assumptions,  it  does  provide  a  good  approximation  of  the  sensitivity  of  the  array  to  errors 
as  a  function  of  signal  direction.  That  is,  the  exact  values  of  the  variance  may  change 
slightly  depending  on  the  error  mechanism  involved,  but  the  shapes  and  trends  (especially 
towards  endfire)  exhibited  by  the  curves  in  Figure  8  would  be  expected  to  stay  the  same. 

Given  that  the  endfire  effects  are  due  solely  to  the  sin“^(-)  transformation  between 
spatial  frequency  and  azimuth  bearing,  and  potentially  could  mask  other  underlying  error 
mechanisms  (besides  noise),  it  is  useful  to  normalize  the  error  results  accordingly.  One 
way  is  to  compute  the  spatial  frequency  errors,  rather  than  the  actual  bearing  errors, 
and  then  convert  these  to  bearing  errors  as  if  the  signal  bearing  had  originally  been  zero 
degrees.  This  is  equivalent  to  using 


^.(t)  =  sin"' 


(35) 


instead  of  equation  (29)  where  t  is  the  spatial  frequency  error.  In  the  c£ise  where  only 
the  true  and  estimated  bearings  are  known,  the  corresponding  spatial  frequencies  can  be 
computed  using  equation  (23).  In  terms  of  accuracy,  the  effect  of  this  normalization  is  to 
artificially  rotate  the  array  so  that  the  array  baseline  is  alway  perpendicular  to  the  signal 
source.  For  any  given  signal  to  noise  ratio,  the  accuracy  is  then  predicted  by  the  values 
for  ^  =  0°  in  Figure  8.  For  obvious  reasons  this  normalization  process  is  called  broadside 
normalization. 


4.3  Wind  Effects 

The  effect  of  the  wind  on  the  Osprey  antenna  array  is  to  perturb  the  antennas  from  their 
known  positions.  Ignoring  all  other  error  mechanisms  the  perturbed  output  at  antenna  n 
will  be, 

j/„  =  (36) 

where  and  Syn  are  the  displacements  of  the  antenna  from  its  ideal  position  in  the  x 
and  y  directions,  respectively.  In  reference  [6],  an  expression  of  the  effects  of  y  direction 
position  errors  (i.e.  ^xn  =  0)  on  bearing  error  variance  is  developed.  In  this  analysis  it  is 
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assumed  that  Syn  has  a  Gaussian  distribution  and  |^y„|  "C  d.  Using  this  expression  and 
converting  it  to  the  equivalent  spatial  frequency  error  variance  (this  is  done  so  that  the 
contribution  in  the  x  direction  can  be  more  readily  included),  the  result  is  given  by 

2  12  f2nd  y 

where  Oy  is  the  position  error  variance  in  the  y  direction  and  is  asstimed  to  be  the  same 
for  all  antennas.  Noting  the  similarity  between  the  effects  of  the  errors  S^n  sin  ^  and 
^y„  cos  <f>  in  equation  (36),  then  one  can  easily  modify  equation  (  37)  to  include  the  effects 
of  0  <  l^xnl  'C  d.  Making  the  same  assumptions  for  Sxn  as  for  Syn  and  also  assuming  Sxn 
and  Syn  to  be  uncorrelated,  the  modified  expression  becomes 

=  n{N^-  1)  (y)  +  (38) 

where  al  is  the  position  error  variance  in  the  y  direction.  In  the  special  case  where 
<Tx  =  <Ty  =  (Txy,  this  simplifies  to, 


2  12  /2wdcrxy\^ 

~  - 1)  V  A  / 


(39) 


The  conversion  of  the  spatial  frequency  error  variance  to  bearing  error  variance  proceeds 
in  the  same  manner  as  the  noise  case  (i.e.  using  equations  (27)-(34)). 

Some  comments  about  the  above  expression  are  in  order  here.  The  antennas  used  for 
the  Osprey  array  have  a  very  small  wind  cross  section  and  are  sxifBciently  rigid  that  for 
all  but  the  most  extreme  weather  conditions  the  values  of  Sxn  and  Syn  are  small  compared 
to  the  antenna  spacing. 

The  antennas  are  also  identical  which  means  that  the  position  error  variances  for  each 
antenna  will  be  identical.  This  also  means  that  the  antennas  would  be  expected  to  bend 
in  approximately  the  same  direction  when  the  wind  blows.  Consequently  the  position 
errors  will  be  correlated  from  sensor  to  sensor.  Since  small  position  translations  of  the 
array  have  no  effect  on  the  DF  results,  then  equations  (38)  and  (39)  can  be  regarded  as 
upper  bounds  on  the  error.  Alternately,  Sxn  and  Syn  (and  correspondingly  <Tx  and  Oy)  can 
be  interpreted  as  the  position  errors  which  result  after  the  effective  position  translation 
of  the  whole  array  due  to  the  wind  has  been  taken  into  account  (i.e.  subtracted  out). 

Over  short  periods  of  time,  the  wind  direction  will  be  constant  with  the  result  that 
the  perturbations  Sxn  and  ^y„  will  be  highly  correlated.  Therefore  equations  (38)  and  (39) 
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axe  assumed  to  apply  for  time  periods  long  enough  that  ^xn  and  Syn  decorrelate. 

The  preceding  analysis  does  not  take  into  account  the  change  in  the  mutual  coupling 
or  polarization  response  of  the  antennas  when  bending  occurs.  These  effects  have  been 
ignored  for  the  following  reasons.  Mutual  coupling  is  caused  by  reradiation  of  the  inci¬ 
dent  electric  field  by  the  antennas  (and  support  structure)  themselves.  Changing  antenna 
position  affects  the  amount  of  reradiated  signal  power  received  at  each  antenna.  Since 
this  is  a  secondary  effect  (i.e.  the  reradiated  power  received  by  one  antenna  from  another 
antenna  will  be  significantly  less  than  from  the  power  received  from  the  transmitter  — 
see  the  following  section  on  mutual  coupling  ),  the  bearing  errors  introduced  by  this  ef¬ 
fect  will  be  correspondingly  less  than  those  discussed  in  the  preceding  paragraphs.  In 
the  case  of  the  polarization  response,  the  amplitude  gain  of  a  dipole  antenna  for  verti¬ 
cally/horizontally  polarized  signals  is  a  function  of  the  cosine/sine  of  the  bending  angle 
[7].  For  the  small  bending  angles  observed  in  practise  (<  1°),  and  even  assuming  the 
incident  vertical  and  horizontal  electric  fields  are  completely  different,  the  bearing  errors 
introduced  are  negligible. 

4.4  Mutual  Coupling 

Mutual  coupling  effects  arise  from  the  fact  that  each  antenna  modifies  the  electric  field 
locally  by  retransmitting  a  portion  of  the  received  signal.  The  amount  of  coupling  be¬ 
tween  antennas  is  therefore  a  function  of  the  ratio  of  the  received  signal  power  to  the 
retransmitted  signal  and  the  distances  between  antennas. 

In  the  case  of  the  Osprey  array,  there  are  several  observations  which  simplify  the  anal¬ 
ysis.  The  first  is  the  small  electrical  radii  of  the  dipole  elements  (<  O.OIA)  which  means 
dipole  current  will  only  be  in  the  vertical  direction.  The  second  is  that  the  directional  or 
beam  pattern  of  a  vertical  dipole  is  omnidirectional  in  azimuth.  The  third  is  that  each 
element  lies  in  the  same  vertical  plane.  All  these  features  ensure  that  an  expression  de¬ 
scribing  mutual  coupling  effects  will  be  dependent  only  on  the  distance  between  antennas 
and  independent  of  signal  bearing  in  both  azimuth  and  elevation. 

Based  on  these  observations,  the  analysis  of  mutual  coupling  can  be  treated  in  terms 
of  an  equivalent  electric  circuit.  For  example,  the  equivalent  circuit  for  antenna  0  in  a 
two  antenna  system  is  shown  in  Figure  9  where  Vo  is  the  voltage  due  the  signal  source 
(transmitter),  Vbi  is  the  voltage  induced  in  antenna  0  via  reradiation  from  antenna  1, 
is  the  self  impedance  of  the  antenna,  Zi^  is  the  load  impedance,  and  Vo  is  the  voltage 
across  the  two  impedances.  The  voltages  in  this  circuit  have  the  relationship 

Vo  =  VoH-Foi  (40) 
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The  voltage  induced  in  antenna  0  by  antenna  1  can  be  related  to  the  current  in  antenna 
1  by  the  expression 

Voi  =  IiZoi  (41) 

where  Zoi  is  defined  as  the  mutual  impedance  between  antennas  0  and  1.  Adopting  a 
similar  equivalent  circuit  definition  for  antenna  1  as  given  for  antenna  0  (shown  in  Figure 
9),  then  the  current  in  antenna  1  can  be  determined  using 


Vi 

Zti  +  Zi^ 


(42) 


Plugging  this  last  result  into  equation  (41)  and  then  back  into  equation  (40)  gives 


Similarly  for  antenna  1, 


(43) 


(44) 
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These  results  are  easily  generalized  to  an  N  element  array  to  get 

(zr+^)  («) 

where  Za  =  Zn  +  Zi^.  Since  Vjt  is  the  voltage  that  would  be  measured  by  an  ideal  sensor 
(i.e.  the  voltage  due  to  the  transmitter  only),  it  is  the  voltage  of  interest  for  DF  purposes. 
In  reality  the  only  sensor  voltages  available  will  be  the  load  voltages  VJ*  (fc  =  0, 1, A^— 1). 
Since  Zi^  and  Zr*  form  a  voltage  divider,  can  be  related  to  14  using 

V,  =  W) 

Equation  (45)  then  becomes 

forA:  =  0,l,2,...,iV-l  (47) 


The  system  of  equations  represented  by  equation  (47)  can  also  be  written  in  matrix 
form  as 

z  =  Qy  (48) 

where  z  represents  the  ideal  antenna  voltage  with  z*  =  V*,  y  represents  the  measured 
antenna  output  voltage  (after  equalization)  with  j/it  =  V/*,  and  Q  is  the  N  x  N  impedance 
ratio  matrix  with  elements 


gki  = 


Zi. 


When  i  =  k  the  above  expression  can  also  be  written  as 

Zvk  +  Zi^ 


gkk 


2i. 


(49) 


(50) 


For  dipole  antennas,  the  mutual  impedance  can  be  represented  as 


Zki^Rki+jXki  (51) 

where  Rki  is  the  resistive  element  and  Xki  is  the  reactive  element.  An  analytical  expression 
for  these  elements  is  given  by  [7],  [8] 
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-4cos‘(y)  [ci(|(V'4rJ.  +  - i))  +  Ci(|( +  i))] 

+  cos(*I)  [Ci(«t(v'r|i  +  +  I))  +  Ci(<t(  V'rL  +  -  I))] 

+  sin(ici)  [si(-c(/tTP  + 1))  -  Si(K(^rl  +  L^-L)) 

-2Si(|(\/4r|i  + +  i))  +  2Si(|(\/4rl,.  +  1^-  I))])  (52) 

—30  / 

(2(2  +  co»(«i))Si(«r„.)- 

-4cos>(!^)  [Si(|(v(4r2,  + -  L))  +  Si(|(V'4r2(  +  +  i))] 

+  cos(Ki)  (si(K(\/r2i  + 1*  + 1))  +  -f  -  X))j 

- sm(/ci)  [ci(K(v'r?,  +  i2  +  i)) _  Ci(«(\Ai+^ -  i)) 

-2Ci(|(y4rJ,  +  +  i))  +  2Ci(|(\/4rJ,  +  -  i))])  (53) 

where  r*,-  =  |A:  —  t|(i  is  the  distance  between  antennas  k  and  t,  k  =  27r/A,  and  Ci(')  and 
Si(-)  are  the  cosine  and  sine  integrals  defined  as 


Ci(a:)  =  r—dv  =  0.5772  +  lnx-^  +  -^--^  +  ... 
.  ^  Joo  V  2!2  4!4  6!6 


(54) 


and 


Si(i)  =  f 
Jo 


^sinv  x°  X 

0  iT  ”  ~  ^  “  3!3 5!5  ~  m " 


(55) 


The  impedance  Zkk  consists  of  the  antenna  self  impedance  Zr*  and  the  load  impedance 
Z/*.  For  maximum  power  transfer  (i.e.  a  matched  system),  the  load  impedance  Z/*  =  Z*^ 
which  results  in  Zkk  =  2real{Zr*}.  The  self  impedance  Z,*  can  be  computed  using 
equations  (51)-(55)  where  rki  is  equated  to  the  radius  of  the  dipole  elements  (0.3  cm  for 
the  dipoles  of  the  Osprey  array). 

Equations  (51)- (55)  predict  identical  mutual  impedance  for  any  pair  elements  provided 
the  spacing  is  the  same.  In  reality,  this  will  not  be  exactly  true  due  to  electrical  and 
physical  diflferences  between  the  antennas,  and  mismatches  introduced  by  cables,  cable 
couplings,  etc.  However  these  equations  should  be  sufiSciently  accurate  to  characterize  the 
resulting  errors  introduced  into  DF  estimation  and  approximately  quantify  these  errors. 
In  this  light.  Figures  10  and  11  show  the  performance  of  the  Osprey  with  mutual  coupling 
effects  included  and  assuming  the  antennas  are  perfectly  matched.  The  frequency  chosen 
for  simulation  was  62.5  MHz. 
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Figure  10:  Estimated  bearing  versus  true  bearings  with  mutual  coupling. 
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Figure  11:  Corresponding  bearing  errors:  (a)  actual  (b)  after  broadside  normalization. 
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Figure  10  plots  the  estimated  bearings  versus  the  true  bearings.  The  most  significant 
errors  occur  towards  the  end  fire  direction  with  the  estimated  bearings  locking  onto  ±90° 
once  the  actual  bearing  reaches  ±77°.  The  lock  on  effect  is  a  result  of  the  mutual  coupling 
which  increases  the  magnitude  of  the  measured  spatial  frequency,  compared  to  the  true 
frequency,  for  signal  directions  near  endfire.  This  causes  a  corresponding  increase  in  the 
measured  bearing  magnitude  so  that  if  the  true  bearing  is  increasing  in  magnitude  the 
measured  bearing  will  reach  ±90°  before  the  true  bearing  does.  Since  a  linear  array 
can  only  measure  bearings  in  the  range  ±90°  to  —90°,  increasing  the  magnitude  of  the 
transmitter  bearing  beyond  this  point  has  no  effect  on  the  measured  bearing,  i.e.  the 
measured  bearing  will  remciin  locked  on  ±90°. 

Figure  11a  shows  the  corresponding  bearing  error.  Figure  11b  shows  the  errors  after 
broadside  normalization  (defined  previously  by  equation  35))  and  illustrates  that  not  all 
the  increased  endfire  error  is  attributable  to  the  conversion  from  spatial  frequency  to 
bearing.  The  RMS  errors  for  Figure  11  were  3.5°  for  (a)  and  0.5°  for  (b). 

A  final,  but  important  observation  about  mutual  coupling  is  that  the  resultant  error 
response  does  not  fluctuate  rapidly  with  azimuth  and  that  it  is  invariant  with  time  (ignor¬ 
ing  equipment  aging  and  minor  wind  effects)  —  important  distinguishing  characteristics 
which  are  useful  when  other  error  effects  are  included  in  the  analysis. 

4.5  Multipath 

VHF  multipath  and  its  effects  on  communications  systems  have  been  extensively  re¬ 
searched  and  reported  in  the  open  literature  (e.g.  IEEE  Transactions  on  Antennas  and 
Propagation).  To  the  author’s  knowledge,  very  little  research  has  been  done  on  VHF 
multipath  and  its  effects  on  direction  finding  antenna  arrays.  DF  research  differs  from 
communications  research  in  the  fact  that  not  only  is  received  signal  power  an  important 
consideration,  but  so  is  the  shape  of  the  incident  signal  wavefront. 

Some  idea  of  the  complexity  of  the  DF  problem  can  be  realized  by  considering  the 
ground  conditions  over  which  the  radio  signal  is  transmitted.  In  the  rural  country  side 
around  DREO,  ground  conditions  are  very  non-uniform,  i.e.,  varying  ground  conductiv¬ 
ities,  numerous  random  obstacles  such  as  rocks,  hills,  trees,  forests,  power  lines,  fences, 
buildings,  etc.,  and  varying  size  and  conductivities  of  these  obstacles.  This  results  in  the 
generation  of  a  multitude  of  secondary  signals  through  diffraction,  reflection,  reradiation, 
refraction  and  scattering  effects,  which  adversely  affects  DF  accuracy. 

There  have  been  efforts  made  at  DREO  to  model  these  effects  with  the  aim  of  better 
imderstanding  the  consequences  for  direction  finding.  Reporting  on  this  research  in  detail 
is  beyond  the  scope  of  this  report,  however,  it  is  worth  describing  some  of  the  more  general 
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results  which  have  relevance  to  the  experimental  measurements  described  in  this  report. 

In  order  to  keep  the  following  analysis  relatively  simple,  a  number  of  assumptions  have 
been  made.  These  assumptions  are: 

1.  Ray  theory  is  applicable 

2.  All  antennas  are  isotropic 

3.  Farfield  conditions  are  assumed 

4.  Antenna  heights  are  insignificant  compared  to  the  transmitter  receiver  range 

5.  The  transmitter  is  the  phase  reference 

The  first  assumption  allows  path  losses  and  corresponding  phase  delays  to  be  calculated 
in  a  relatively  straight  forward  manner.  However,  it  does  not  provide  any  indication  of 
the  scattering  that  occurs  when  the  signal  wavefront  encounters  an  obstacle.  The  second, 
third,  and  fourth  assumptions  simplify  the  mathematics  and  only  degrade  the  results 
when  the  effect  of  multipath  sources  within  a  few  wavelengths  of  either  the  transmitting 
or  receiving  antennas  are  considered.  The  last  assumption  provides  a  common  reference 
for  the  phase  calculations. 

In  the  following  three  sections,  the  multipath  model  is  developed  beginning  with 
freespace  propagation,  followed  by  the  inclusion  of  ground  reflections,  and  ending  with 
the  development  of  idealized  multipath  sources  and  their  effect  on  propagation. 

4.5.1  Freespace  model 

Proceeding  from  the  assumptions  just  described,  the  received  signal  can  be  determined  for 
a  given  transmitted  signal  based  on  power  considerations.  For  an  isotropic  transmitter 
source  in  freespace,  the  radio  wave  expands  in  a  spherical  wave  with  the  transmitter 
source  at  the  center.  The  power  extracted  by  a  receiving  antenna  at  a  distance  r  from 
a  transmitting  antenna  will  be  in  direct  proportion  to  the  area  of  the  spherical  wave  ft 
intercepted  by  the  antenna  compared  to  the  total  area  47rr*  of  the  spherical  wave  —  see 
Figure  12.  If  the  collecting  area  of  the  antenna  is  circular,  the  received  power  will  be 
given  by. 


where  Prew  is  the  received  signal  power,  Pn  is  the  power  radiated  by  the  transmitter, 
and  ft  is  the  effective  collecting  area  of  the  receiving  antenna.  For  an  isotropic  receiving 
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antenna  ^  ^  [7],  and  equation  (56)  becomes, 


For  r  X  this  simplifies  to 


(57) 


(58) 


It  is  often  useful  to  express  the  path  loss  as  a  relative  quantity.  Accordingly,  the 
freespace  path  loss  Ijrtc,  is  defined  here  as  the  ratio  of  the  received  signal  power  to  the 
transmitted  power,  or 


/ 


'free 


=  ( 


Airr  ) 


(59) 


The  complex  amplitude  Snvr  of  the  received  signal  can  also  be  related  to  that  of  the 
amplitude  of  the  transmitted  signal  S(j  (which  is  real  valued  since  the  transmitter  is  the 
phase  reference)  using  the  fact  that  |s|  a  Pa  (where  s  and  P  represent  signal  amplitude 
and  power  respectively)  and  the  phase  delay  is  a  function  of  the  path  length  r.  This  leads 
to  the  result 

/  \  \ 

(60) 


-(sfc) 


;  ^  A 
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4.5.2  Plane  Earth  Model 


For  signal  paths  over  land,  the  ground  plays  a  prominent  role  in  propagation  loss.  Begin¬ 
ning  with  the  plane  earth  model,  as  shown  in  Figure  13,  the  received  signal  is  the  sum 
of  the  direct  signal  and  the  ground  reflected  signal.  The  direct  signal  can  be  calculated 
using  the  freespace  equations  to  give 

(61) 

where  st^  is  the  amplitude  of  the  signal  measured  at  the  transmitter.  The  ground  reflected 
signal  can  be  calculated  in  a  similar  fashion  to  get 

where  Ry  is  the  ground  reflection  coefficient.  In  terms  of  the  direct  signal  Sd,>  this  becomes 

Sref  =  SdirR,  (63) 

Vri-i-r2/ 


Figure  13:  Plane  Earth  propagation  model. 


Assuming  ri,r2  hi,h2  then  the  first  br2w:keted  term  in  equation  (63)  can  be  ap¬ 
proximated  as 

rp 

ri  +  r2 

and  further  assuming  rX  ^  h\h2  then  the  second  bracketed  term  (which  is  part  of  the 
exponential  term)  in  equation  (63)  can  be  approximated  by 


ri  +  r2  -  rp  = 


2hih2 

r 


(65) 
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are  valid,  and  the  expression  for  the  reflected  signal  s^e/  can  be  simplified  to 


SreJ  ~  ^dirRv^  ^ 


^dirRv 


(.  Anhih2\ 

j 


The  combined  signal  Srcur  =  ■Sav  +  Srej  is  given  by, 


^rcvT  ““  j Rv 


4whih2\ 


(66) 

(67) 


(68) 


Figure  14:  Reflection  coefficient  versus  grazing  angle. 

The  reflection  coefficient  is  a  fimction  of  several  factors.  For  vertically  polarized  signals 
the  appropriate  expression  is  given  by, 

Cc  sin  tf}  —  y/tc  —  cos^  rp 


Rr,= 


£c  sin  Ip  +  y/tc  —  cos^  rp 


(69) 
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where  ^  is  the  grazing  angle  which  can  be  determined  using 


V,  =  sin-  (70) 

and  Cc  is  the  complex  dielectric  constant  which  can  be  defined  in  terms  of  the  dielectric 
constant  c,  and  conductivity  o  of  the  reflection  surface  as  [9] 


Cc  =  Cr— iflOffA  (71) 

Figure  14  plots  the  real  and  imaginary  values  of  Ry  for  dry  earth  and  wet  conditions. 

Considering  the  assumptions  made  for  r,  hi,  and  hi  the  reflection  coefficient  can  be 
approximated  by 

2c 

Rv  »  sin  ^  -  1  (72) 

V^c  —  1 

Substituting  the  right  hand  side  of  equation  (70)  for  %})  and  using  the  resultant  expression 
to  replace  Ry  in  equation  (68)  then 


,,  Sdir  /2cc(hi  +  hj)  ,  .47rhih2\ 
r  V  VCc^  A  J 

In  terms  of  the  transmitted  signal  (equation  (61)) 


(73) 


“  Ta"  +  ^2)  +  7^1  fta)  (74) 

where  the  fact  that  tq  r  has  been  used  here.  For  the  frequencies  and  antenna  heights 
used  during  the  field  trials,  the  above  expression  for  is  accurate  to  within  0.2  dB 
amplitude  and  0.3°  phase  at  a  range  of  1  km,  and  4  dB  amplitude  and  6°  phase  at  a  range 
of  50  m. 

The  corresponding  path  loss,  Icarth  can  be  calculated  using 

=  (75) 

Approximating  Srevr  by  equation  (74),  the  path  loss  expression  becomes 


_i7P 

earth  —  ^ 


(76) 
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where 

^  =  (77) 

4.5.3  Multipath  Model 

In  examining  Figure  13,  the  signal  is  reflected  off  the  ground  at  a  single  point.  In  reality, 
the  entire  ground  contributes  to  the  reflected  signal  (i.e.  each  point  on  the  siufface  can 
be  considered  as  a  point  source  reradiator).  It  is  the  assumption  of  a  smooth  uniform 
infinite  earth  that  leads  to  a  reflected  wave  which  is  a  mirror  image  of  the  incident  wave, 
and  therefore  allows  the  use  of  ray  theory. 

In  most  practical  situations,  the  ground  will  have  features  which  cause  it  to  depart 
from  uniform  smooth  earth  conditions  (e.g.  terrain  features  and/or  areas  with  different 
values  of  Cc).  The  result  will  be  a  change  in  the  signal  which  is  reradiated  from  the 
ground  surface  of  each  feature.  The  approach  used  here  is  to  consider  the  reradiated  signal 
from  the  feature  as  the  sum  of  a  smooth  earth  component  plus  a  multipath  component. 
The  cumulative  effects  of  the  smooth  earth  signals  from  the  entire  earth  surface  can  be 
modelled  using  the  plane  earth  model  and  are  therefore  ignored  in  the  rest  of  this  analysis. 
It  only  remains  necessary  to  develop  a  model  for  the  multipath  component. 

For  simplicity  the  features  which  cause  a  departure  from  uniform  smooth  earth  con¬ 
ditions  are  called  multipath  sources.  Since  multipath  sources  are  not  infinite  planes,  the 
reradiation  of  the  signal  (multipath  component  only)  from  a  multipath  source  will  re¬ 
sult  in  some  scattering  (spreading)  of  the  reradiated  wavefront  compared  to  the  incident 
wave.  Quantifying  this  effect  is  beyond  the  scope  of  this  report  so  that  only  the  two  ex¬ 
treme  scattering  conditions  are  considered:  no  scattering  and  isotropic  scattering  (i.e.  the 
multipath  feature  acts  like  an  isotropic  source).  For  convenience,  a  multipath  source  for 
which  no  scattering  occurs  will  be  called  a  perfect  reflector,  and  a  source  which  produces 
isotropic  scattering  will  be  called  an  isotropic  scatterer.  The  bearing  of  the  reradiated 
signal  may  also  be  different  from  that  of  the  incoming  signal  due  to  the  underlying  mech¬ 
anisms.  Although  this  fact  is  not  overlooked  in  the  following  analysis,  no  attempt  is  made 
to  quantify  this  effect. 

Figiire  15  illustrates  the  case  for  a  single  perfect  reflector  as  viewed  from  above  and 
the  side.  This  case  is  nearly  identical  to  the  direct  signal  case,  except  that  the  “direct” 
signal  travels  from  the  transmitter  to  the  the  multipath  source  (ro)  and  on  to  the  receiver 
(rj).  Viewed  from  above,  the  signal  path  will  not  be  a  straight  line  unless  the  transmitter, 
multipath  source,  and  receiver  are  all  aligned.  Likewise,  the  groimd  reflected  signal  travels 
from  the  transmitter  to  the  ground  (rj),  from  the  ground  to  the  multipath  source  (r2a), 
and  from  the  multipath  source  to  the  receiver  (rjb).  Note  that  in  this  example,  the 
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I  I 


(b)  Side  view 

Figure  15:  Plane  Earth  propagation  model  for  a  perfect  reflector. 

The  multipath  source  will  delay  and  attenuate  the  signals,  requiring  the  introduction 
of  a  complex  radiation  coefficient  Rm  where  |ilnj|  <1.  It  is  assumed  here  that  is 
identical  for  both  signal  components. 

Given  the  similarity  between  this  case  and  the  plane  earth  model,  and  using  the  same 
assumptions,  equations  (74)  and  (76)  are  easily  modified  to  get 


(ra  +  ri)2 


and 


im%H  — 


\lRm? 


(78) 


(79) 


{va  +  n)* 

It  is  also  useful  to  compare  the  path  loss  in  this  case  to  the  path  loss  leartk  undergone 
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by  the  direct  signal.  Defining  the  ratio  lexers  =  imuhlleartk  as  the  excess  path  loss,  then 

=  (80) 
\ra  +  nJ 

Figure  16  plots  the  excess  path  loss  as  a  function  of  position  for  the  transmitter  and 
receiver  separated  by  5  km.  It  is  assumed  that  for  each  position  on  the  plot  the  geometry 
is  favourable  for  the  generation  of  a  reradiated  signal  directed  towards  the  receiver.  Since 
J2m  simply  scales  the  results,  a  value  of  \Rm\  =  1  is  used.  Not  surprisingly  there  is  no 
excess  path  loss  for  multipath  sources  located  along  the  direct  signal  path  with  path  loss 
increasing  as  the  position  is  moved  from  the  direct  signal  path. 

Figure  17  illustrates  the  case  for  an  isotropic  scatterer.  Since  the  multipath  feature 
acts  like  an  isotropic  source,  the  path  losses  from  the  transmitter  to  the  multipath  source 
and  from  the  multipath  source  to  the  receiver  can  be  treated  independently.  Hence  the 
received  signal  is  given  by. 


where 


S..U  = 

(81) 

(82) 

(83) 

and  hm  is  the  height  of  the  signal  path  above  the  ground  at  the  multipath  source  (shown 
Figure  17).  Simplifying  the  expression  for  Sm*h  leads  to 


^rnnh 


(84) 


The  corresponding  path  loss  equation  is  given  by 

I-Rm7i72p 


Im»/<  — 

The  excess  path  loss  in  this  case  is  given  by, 

iexeesi 


{ranY 


(85) 
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(86) 


Figure  18  plots  the  excess  path  loss  for  the  same  transmitter-receiver  geometry  and 
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(b)  Side  view 

Figure  17:  Plane  Earth  propagation  model  for  an  isotropic  scatterer. 
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assumptions  as  used  in  Figure  16.  The  ground  conductivity  was  also  assumed  to  be 
<7  =  3  X  10“^ JWho/m,  and  the  dielectric  constant  =  3  (dry  earth  conditions).  For 
simplicity,  the  center  of  the  scattering  source  was  assumed  to  lie  on  the  line  joining  the 
transmitting  and  receiving  antennas  when  the  horizontal  profile  is  viewed  as  shown  in 
Figure  17.  Accordingly, 

,  2hih2 

ir+hi 

However,  in  general,  km  will  depend  on  the  size  and  shape  of  the  feature  so  that  the 
above  expression  will  not  usually  be  true.  For  positions  corresponding  to  the  transmitter 
and  receiver,  it  was  assumed  that  there  was  no  excess  path  loss.  Inspection  of  Figure 
18  reveals  that  in  the  scattering  case  the  excess  path  loss  increases  as  a  function  of  the 
distance  from  both  the  transmitter  and  receiver,  not  the  direct  signal  path  as  for  the 
reflector  case  shown  in  16.  The  loss  also  increases  at  a  faster  rate  for  the  scattering  case 
than  for  the  reflector  case. 

Before  drawing  too  many  conclusions  from  a  comparison  of  these  results,  it  is  worth 
considering  another  factor  which  is  important:  the  direction  of  the  multipath  source 
relative  to  the  transmitter.  For  example,  if  the  multipath  source  is  aligned  with  the 
transmitter  and  is  also  located  in  the  farfield  of  the  receiving  array,  the  overall  effect 
will  be  to  change  the  power  of  the  received  signal,  but  not  the  relative  signal  phases  and 
amplitudes  measured  at  the  array.  Since  the  bearing  estimates  depend  only  on  the  relative 
amplitude  and  phase  information,  no  errors  are  introduced.  If  the  multipath  source  is  in 
a  different  direction  than  the  transmitter,  then  this  will  no  longer  be  the  case. 

In  an  attempt  to  quantify  this  effect,  and  ignoring  all  other  errors  except  multipath, 
a  single  snapshot  of  equalized  data  can  be  defined  as 

y  =  /5s  +  m  (88) 

where  /5  is  the  complex  signal  amplitude,  s  is  the  signal  vector  defined  by  equation  (12), 
and  m  is  the  IV  x  1  vector  representing  the  contribution  of  multipath.  The  component  of 
multipath  that  gives  rise  to  error  in  the  DF  estimates  is  that  which  is  orthogonal  to  the 
signal  vector  s.  Therefore  it  is  useful  to  redefine  equation  (88)  as, 

y  =  Qs  +  mj.  (89) 

where  a  is  the  apparent  signal  amplitude  given  by 

Q  =  ^  +  B^m  (90) 
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Multipath  Attenuation  (dB) 


and  mx  is  the  orthogonaJ  component  of  the  multipath  such  that  s^mx  =  0  and  is  given 

by, 

mx  =  m  —  ss^m  (91) 

Defining  the  loss  function  as  the  ratio  of  the  apparent  multipath  power  to  the  actual 
multipath  power,  and  calling  this  the  directional  loss  function,  then 


= 


Imp  —  Im^sp 


=  1-|PP 


For  a  single  multipath  source  in  the  farfield  of  the  receiving  array  m  will  have  the  form 
given  in  equation  (12)  except  the  transmitter  bearing  <f>  is  replaced  by  the  multipath 
source  bearing  (p.  Under  these  conditions 

^  N{1  — 

which  simplifies  to  p  =  1  for  (^  =  ^. 

The  function  is  plotted  (in  dB)  in  Figure  19  for  the  same  transmitter-receiver 
geometry  used  in  Figure  19  and  assuming  the  direct  signal  bearing  is  ^  =  0°.  As  would 
be  expected,  multipath  sources  approximately  in  line  with  the  transmitter  and  receiver 
have  far  less  effect  than  sources  which  are  further  off  the  transmitter-receiver  line. 

It  is  useful  to  observe  the  combined  effect  of  and  lexce»$  for  both  types  of  multipath 
considered.  Accordingly,  Figures  20  and  21  are  repeats  of  Figures  16  and  18  with  the 
directional  effects  included.  The  main  difference  between  Figures  20  and  21  and  the 
originals  is  the  suppression  of  multipath  originating  near  the  transmitter,  although  this 
is  far  more  evident  in  the  case  of  the  scattering  somce  of  multipath. 

Based  on  these  figures  some  further  comments  are  in  order.  For  example,  perfect 
reflectors  are  geometry  dependent  —  they  illuminate  the  receiving  array  for  a  very  limited 
subset  of  transmitter  bearings  and  ranges.  The  result  is  that  the  DF  bearing  error  for 
one  transmitter  bearing  will  be  completely  unrelated  to  the  error  occurring  at  another 


(94) 
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Figure  20:  Combined  effect  of  excess  path  loss  and  directional  loss  for  a  perfect  reflector. 
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bearing.  This  effect  is  illustra,ted  in  Figure  22  where  the  DF  error  response  for  the  Osprey 
system  (after  broadside  normalization  as  defined  in  equation  (35))  has  been  simulated  for 
a  single  signal  and  three  multipath  sources  at  equal  ranges.  Each  source  was  assumed 
to  illuminate  the  receiving  array  only  when  —  <j>  =  10°,  where  m  is  the  multipath 
source  designation  (1,  2,  or  3).  The  three  plots  show  the  effect  of  increasing  the  distance 
between  the  receiver  and  the  multipath  sources  from  10  meters  to  1000  meters.  Note  that 
the  results  at  close  range  are  unrealistic  in  the  sense  that  a  perfect  reflector  will  have 
infinite  size  which  implies  near  field  conditions,  i.e.  simple  ray  tracing  will  not  be  valid. 
However  the  results  do  provide  a  reasonable  indication  of  the  effects  of  both  range  and 
bearing. 


Figure  22:  Simulated  DF  error  response  for  3  multipath  reflectors  at  bearings  of  -60,  40, 
and  90  degrees,  and  ranges  of  (a)  10  m  (b)  100  m,  and  (c)  1000  m. 


Examining  the  plots  in  Figure  22,  the  effects  of  the  individual  multipath  sources  are 
clearly  identifiable  (at  <f>  =  —70,  30,  and  80  degrees)  and  uncorrelated.  Range  does  not 
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Figure  23:  Simulated  DF  error  response  for  3  multipath  scatterers  at  bearings  of  -60,  40, 
and  90  degrees,  and  ranges  of  (a)  10  m  (b)  100  m,  and  (c)  1000  m. 
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appear  to  have  any  significant  effect.  Extrapolating  these  results  to  a  greater  number 
of  multipath  sources  would  result  in  a  very  noise-like  error  pattern  with  no  correlation 
between  the  error  at  any  one  bearing  compared  to  any  other  bearing.  The  implication 
is  that  if  a  large  number  of  reflector  like  sources  dominate  the  multipath  environment, 
then  calibration  of  the  array  response  will  be  very  difficult,  i.e.,  incremental  measurements 
will  be  required  over  the  complete  spectrum  of  transmitter  azimuth  angles  and  ranges  of 
interest  to  achieve  any  significant  accuracy  improvement  compared  to  equalization. 

Isotropic  scattering  sources  are  not  dependent  on  the  bearing  of  the  incident  signal 
but  are  affected  by  the  path  lengths  involved.  The  strongest  multipath  components  will 
tend  to  originate  near  the  receiving  antenna.  If  the  transmitter  is  moved,  the  multipath 
environment  will  not  significantly  change.  As  a  result,  characterizing  the  site  around 
the  receiving  array  may  be  sufficient  in  order  to  significantly  improve  DF  performance. 
Examples  of  this  are  shown  in  Figure  23  which  is  identical  to  Figure  22  except  that 
the  multipath  sources  completely  scattered  the  signal.  In  this  example  it  is  clear  that 
multipath  effects  are  localized  to  the  area  near  the  receiving  array  and  diminish  rapidly 
with  range.  Additionally,  since  the  error  pattern  is  affected  at  all  azimuth  angles,  this 
implies  a  high  degree  of  correlation  between  measurements  taken  at  different  azimuth 
angles.  A  final  observation  is  that  the  error  pattern  becomes  more  noise-like  as  the  range 
increases,  thus  providing  one  way  of  distinguishing  multipath  effects  from  mutual  coupling 
effects. 


5.0  RESULTS 

The  DF  results  for  all  four  data  sets  are  shown  plotted  in  Figures  24-27.  In  each  figure 
the  top  plot  shows  the  bearing  measurements  ^  (dots)  and  map  bearing  ^  (solid  line) 
in  degrees  as  a  function  of  the  measurement  number.  The  measured  results  have  been 
corrected  for  the  directional  ambiguities  inherent  in  the  antenna  array  (i.e.  a  linear  array 
cannot  distinguish  between  the  two  angles  ^  and  180°  —  if).  The  middle  plot  shows  the 
bearing  error  ^  ^  (in  degrees).  The  lowest  plot  shows  the  bearing  error  (in  degrees) 

after  broadside  normalization  (defined  in  equation  (35)). 

5.1  Error  Analysis 

The  RMS  error  for  each  of  the  data  sets  displayed  in  Figures  24-27  is  given  in  Table 
1  for  both  the  unnormalized  and  broadside  normalized  cases.  Inspection  of  the  figures 
shows  that  the  unnormalized  results  are  dominated  by  the  errors  that  occur  for  bearings 
towards  endfire  (±90°),  especially  for  the  data  sets  collected  in  March.  By  contrast,  the 
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Figure  24:  Field  trial  results  for  November  5,  1992 
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Figure  25:  Field  trial  results  for  November  6,  1992 
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Measurement  number 

Fig;ire  26:  Field  trial  results  for  March  24,  1993 
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Figure  27:  Field  trial  results  for  March  25,  1993 
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normalized  results  exhibit  far  less  dependency  on  DF  antenna  array  geometry  making  this 
representation  of  the  errors  more  suitable  for  analysis.  Therefore  the  remaining  discussion 
of  the  results  deals  with  the  broadside  normalized  results  only. 

In  the  next  three  sections,  the  error  results  are  broken  down  into  temporal  effects, 
spatial  effects  excluding  multipath,  and  multipath  effects.  Figures  28-31  show  the  cor¬ 
responding  range  and  signal  to  noise  ratios  for  the  field  trial  measurements  which  are 
relevant  to  the  following  analysis.  The  broadside  normalized  error  (which  was  shown  in 
Figures  24-27)  is  also  shown  for  reference. 

The  bearing  error  for  measurement  109  on  March  24  is  discussed  in  more  detail  under 
multipath  effects. 

5.1.1  Temporal  Effects 

Temporal  effects  include  noise,  wind  effects,  and  receiver  drift.  Noise  and  wind  effects 
can  be  quantified  based  on  the  analysis  presented  in  Section  4.  Receiver  drift  can  be 
quantified  based  on  the  measurements  as  will  be  discussed  shortly. 

From  Figures  28c-31c  the  signal  to  noise  power  ratio  was  normally  than  40  dB  in  most 
of  the  DF  measurements  during  the  field  trials.  The  associated  RMS  bearing  errors  given 
by  equation  (39),  which  ignores  the  wind  effects  and  receiver  drifts  was  therefore  less 
than  0.1°  and  can  be  considered  insignificant  compared  to  the  RMS  error  values  listed  in 
Table  1.  In  the  few  cases  where  the  signal  to  noise  power  ratio  was  less  than  40  dB  (e.g. 
measurements  9,  10,  79,  and  80  in  Figure  30c),  these  measurements  can  be  ignored. 

To  minimize  wind  induced  bearing  errors,  calm  days  (i.e.  windspeed  less  than  10 
km/h)  were  chosen  for  the  field  trials.  On  these  days  the  position  error  of  the  antennas 
in  the  DF  array  was  less  than  1  cm.  This  translates  to  a  maximum  RMS  bearing  error  of 
0.1°  which  indicates  wind  effects  can  be  ignored. 

To  validate  the  noise  and  wind  results,  27  repeated  DF  measurements  spanning  a 
period  of  3.5  minutes  were  taken  for  the  last  transmitter  position  on  March  25.  The 


Table  1:  RMS  Error  Summary  of  Field  Trial  Results 


Data  Set 

RMS  Error 

Normalized  RMS  Error 

Nov  5 

0.64 

0.51 

Nov  6 

0.72 

0.54 

Mar  24* 

4.44 

1.06 

Mair  25 

4.84 

1.15 

*  Bearing  error  for  measurement  109  excluded 
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(a) 


(b) 


(c) 

Figure  29:  Field  trial  parameters  for  November  6,  1992:  (a)  broadside  normalized  error, 
(b)  range,  and  (c)  signal  to  noise  ratio. 
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Figure  30:  Field  trial  parameters  for  March  24,  1993:  (a)  broadside  normalized  error,  (b) 
range,  and  (c)  signal  to  noise  ratio. 
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Figure  31:  Field  trial  parameters  for  March  25,  1993:  (a)  broadside  normalized  error,  (b) 
range,  and  (c)  signal  to  noise  ratio. 
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standard  deviation  in  these  results  was  0.08°  supporting  the  conclusion  that,  for  high 
signal  to  noise  power  ratios  (>  40  dB)  and  low  winds  (<  10  km/h),  the  noise  and  wind 
induced  errors  were  insignificant. 

Tests  of  the  Osprey  receivers  in  the  lab  found  them  to  be  very  stable  when  properly 
warmed  up  with  phase  drift  rates  of  less  than  0.1°/12  hours.  Although  this  implies 
insignificant  drift  in  the  DF  results,  given  that  the  drift  tests  were  done  under  constant 
temperature  conditions,  the  lab  results  don’t  guaraintee  stable  operation  in  the  field.  To 
test  whether  any  significant  drift  occurred  during  the  measurements,  the  complex  gain 
coeflBcients  were  estimated  from  data  collected  for  the  same  positions  (or  nearly  the  same 
positions)  but  at  different  times,  and  the  results  compared.  There  are  two  groups  of  3 
data  subsets,  for  a  total  of  6  subsets,  in  which  positions  were' repeated.  These  subsets 
are  siunmarized  in  Table  2  with  subsets  A,  B,  and  C  representing  one  group  of  positions, 
and  subsets  D,  E,  and  F  representing  the  second  group  of  positions.  The  time  column  in 
Table  2  refers  to  the  time  of  the  middle  measurement  of  the  subset  (i.e.  the  mean  time). 


Table  2:  Data  Subsets  for  Drift  Estimation 


Subset 

Date 

Time 

Measurement  No. 

Comment 

A 

15:34:55 

all 

see  Figure  3 

B 

12:16:50 

1-81 

same  as  A  (see  also  Figure  4) 

C 

Nov  6 

13:26:53 

115-195 

same  as  A  except  in  reverse  order 

D 

Mar  24 

13:31:55 

13-56 

(xe-Wj  yN-s)  positions  starting  at 
(50,  -472)  to  (-495,  -143)  in  Figure  5 

E 

Mar  24 

14:16:49 

57-102 

same  as  D  except  in  reverse  order 

F 

Mar  25 

10:39:33 

52-197 

same  as  D  (see  also  Figure  4) 

The  estimated  gains  and  phases  for  each  of  the  subsets  are  listed  in  Tables  3  and  4. 
Inspection  of  these  results  would  suggest  that  during  same  day  measurments,  phase  and 
gain  drifts  did  occur,  with  a  maximum  phase  drift  of  up  to  8°  (Table  3,  channel  2)  and 
gain  drift  of  up  13%  (e.g.  Table  3,  channel  4).  These  values  are  substantially  greater  than 
what  would  have  been  predicted  from  the  lab  tests.  However,  since  the  phase  and  gain 
estimates  have  been  constrained  according  to  equations  (8)  and  (9),  the  results  in  Tables 
3  and  4  do  not  necessarily  reflect  the  true  drifts.  For  example,  since  all  the  channel  phases 
are  referenced  to  channel  0,  a  change  in  the  phase  offset  of  channel  0  will  appear  as  a 
constant  phase  change  in  all  the  other  channels.  For  this  reason  interpreting  the  results 
in  this  form  may  be  misleading. 

Additionally,  to  observe  drift  effects  in  more  detail,  DF  measurements  were  taken  over 
periods  of  several  days  with  the  transmitter  at  close  range  and  in  a  fixed  position.  This 
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led  to  the  identification  of  equipment  problems  (i.e.  shielding  problems  which  had  no 
effect  on  the  resiJts  of  the  other  experiments)  which  were  subsequently  fixed.  Although  it 
was  not  possible  to  run  another  properly  controlled  drift  test  for  this  report,  preliminary 
testing  indicates  that  the  environment  may  cause  daily  phase  variations  on  the  order  of 
a  few  degrees.  However,  until  further  tests  can  be  performed  to  confirm  this,  it  has  been 
assumed  that  the  equipment  was  the  sole  cause  of  the  drift  effects. 


Table  3:  Channel  Phase  Drift  during  Field  Trials  in  November 


Channel 

Oa 

6b 

6c 

9A 

9b 

9c 

0 

0.0 

0.0 

0.0 

Era 

0.88 

1 

98.2 

104.7 

104.7 

IQ 

1.02 

2 

44.1 

48.6 

49.3 

IQ 

0.95 

1.03 

3 

101.2 

102.1 

110.1 

IrIiI 

0.94 

0.97 

4 

-167.2 

-161.3 

-161.0 

0.95 

1.17 

1.04 

5 

158.6 

171.1 

170.9 

1.17 

1.30 

1.25 

6 

-88.0 

-75.3 

-75.2 

0.57 

0.44 

0.45 

7 

118.4 

122.0 

118.3 

1.00 

1.12 

1.16 

Table  4:  Channel  Phase  and  Gain  Drift  during  Field  Trials  in  March 


Channel 

6d 

6e 

Bp 

9d 

9e 

9F 

0.0 

0.0 

0.0 

Era 

Era 

-44.6 

-42.7 

-46.4 

Hi 

-35.1 

-35.2 

-39.1 

1.05 

IQ 

0.94 

-60.2 

-59.3 

-65.8 

0.71 

ImI 

0.72 

-72.6 

-68.6 

-75.9 

1.16 

1.20 

1.25 

5 

-23.6 

-21.1 

-27.0 

1.36 

1.36 

1.28 

6 

-5.2 

-3.3 

-9.4 

0.73 

0.71 

0.77 

7 

-14.5 

-10.8 

-12.6 

0.89 

0.87 

0.90 

To  provide  a  better  idea  of  the  effects  of  drift,  it  is  simpler  to  observe  the  error  that  is 
introduced  into  the  DF  results.  Mathematically  the  effects  of  drift  on  a  single  snapshot 
of  data  z  can  be  represented  by 

z  =  Dy  (95) 

where  D  is  an  A"  x  diagonal  drift  matrix  and  y  represents  the  calibrated  data.  The 
diagonal  elements  of  the  drift  matrix,  dQ,di,...,dN,  represent  the  change  in  the  complex 
gain  coefficients  of  the  corresponding  channels.  The  effect  on  DF  estimation  can  therefore 
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Table  5:  Bearing  Errors  (in  degrees)  due  to  Phase  and  Gain  Drift 


Subset 

Bearing  Error 

A 

-0.38 

B 

0.00 

C 

o 

1 

D 

-0.17 

E 

0.00 

F 

-0.16 

be  observed  by  comparing  the  bearing  estimates  determined  using  y  and  z.  Additionally 
if  y  represents  a  single  signal  with  a  bearing  of  0°,  then  the  value  of  the  bearing  estimate 
for  z  will  represent  the  bearing  error  due  to  drift. 

Using  subsets  B  and  E  as  references  for  each  group,  the  diagonal  elements  of  the  drift 
matrix  D  were  calculated  for  each  subset  based  on  the  differences  between  it  and  the 
appropriate  reference  subset.  For  example,  for  subset  A 

-  i^c^iSAk-eBk)  (95) 

9Bk 

where  k  represents  the  channel  number.  Computing  z  for  each  subset  assuming  y  repre¬ 
sented  a  signal  with  a  bearing  of  0°  and  no  noise,  the  estimated  bearings  were  determined 
using  MUSIC  and  tabulated  in  Table  5.  Although  these  results  are  limited,  they  indicate 
that  for  same  day  measurements  the  effect  of  receiver  drift  on  the  bearing  estimates  was 
less  than  0.2°.  This  provides  some  confidence  that  receiver  drift,  aJthough  significant,  was 
not  a  major  contributing  error  during  the  field  trials.  For  consecutive  day  results  the  drift 
rate  was  less  than  0.4°/day,  thus  providing  a  further  indication  that  drift  was  not  a  major 
factor  during  measurements  performed  on  the  same  day. 

5.1.2  Spatial  Effects 

Two  sources  of  spatial  error  are  considered  here,  namely,  transmitter  position  errors,  and 
mutual  coupling.  Both  effects  can  be  quantified  based  on  the  analysis  presented  in  Section 
4. 

Uncertainty  in  the  transmitter  position  degrades  the  accuracy  in  the  determination  of 
the  transmitter  bearings  from  the  map,  especially  as  the  transmitter  gets  closer  to  the  DF 
array.  During  the  November  trials  ranges  varied  from  5.3  to  7.0  km.  The  associated  RMS 
map  bearing  errors  (calculated  using  equation  (21))  ranged  from  0.03°  to  0.02°  RMS, 
which  is  significantly  less  than  the  overall  RMS  bearing  error  values  listed  in  Table  1. 
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Table  6:  Simulated  Bearing  Errors  due  to  Mutual  Coupling 


Data  Set 

Without  Equalization 

With  Equalization 

Nov  5 

0.37 

0.19 

Nov  6 

0.39 

0.23 

Mar  24 

0.73 

0.68 

Mar  25 

0.67 

0.67 

Therefore  for  the  November  trials,  errors  in  the  map  bearing  estimates  were  assumed  to 
be  negligible.  During  the  March  trials,  ranges  were  much  shorter  and  varied  from  22.4  m 
to  1.4  km.  The  associated  RMS  map  bearing  errors  were  7.3°  to  0.15°.  The  larger  RMS 
errors  are  clearly  significant.  However,  by  restricting  analysis  to  measurements  collected 
for  ranges  of  400  m  or  more,  the  resultant  error  in  the  map  bearing  estimates  (<  0.2°) 
can  be  ignored. 

In  the  case  of  mutual  coupling  the  predicted  error  based  on  assuming  perfect  matching, 
and  plotted  in  Figure  11,  was  0.5°  RMS.  This  clearly  is  significant.  To  determine  whether 
the  theoretical  effects  predicted  were  actually  observed  during  the  field  trials,  the  curve 
from  Figure  10  has  been  replotted  in  Figure  32  along  with  the  DF  results  from  March 
25.  In  this  case  the  directional  ambiguities  inherent  to  the  linear  antenna  array  were  not 
corrected  leaving  the  bearings  restricted  to  the  range  —90°  to  +90°.  Inspection  of  the 
results  shows  that  there  is  a  very  good  correlation  between  the  actual  and  theoretical 
results. 

Given  these  good  results,  the  mutual  coupling  simulation  wets  run  for  the  same  trans¬ 
mitter  directions  and  frequencies  used  in  the  actual  trials.  The  RMS  error  results  are 
tabulated  in  Table  6  in  the  middle  column.  Not  smprisingly,  the  March  results,  which 
had  a  greater  number  of  measurements  of  signals  in  the  endfire  direction,  show  greater 
errors  than  the  November  results.  The  stochastic  calibration  method  was  also  applied  to 
the  simulated  data  to  determine  the  combined  effect  of  mutual  coupling  and  equalization. 
The  results  are  also  tabulated  in  Table  6  with  the  November  data  showing  the  greatest 
benefit  from  equalization.  Based  on  these  results,  it  seems  prudent  to  treat  measurements 
of  signals  in  the  endfire  direction  with  some  caution. 

5.1.3  Multipath  Effects 

To  minimize  the  effects  of  the  error  sources  discussed  in  the  previous  two  sections,  the 
measurement  data  was  reprocessed.  Measurements  were  rejected  if  the  signal  to  noise  ratio 
was  <  40  dB,  or  the  transmitter  range  was  <  400  m.  Imposing  this  criteria  resulted  in 
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Figure  32:  Theoretical  mutual  coupling  effects  (solid  line)  compared  to  results  from  March 
25,  1993  (dots) 
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Table  7:  RMS  Error  Summary  of  Field  Trial  Results  with  Mutueil  Coupling  Effects  Re¬ 
moved 


Data  Set 

RMS  Error 

Nov  5 

0.47 

Nov  6 

0.47 

Mar  24* 

0.84 

Mar  25 

0.80 

*  Bearing  error  for  measurement  109  excluded 

2%  of  the  November  measurements  being  discarded  and  16%  of  the  March  measurements 
being  discarded.  The  error  results  are  shown  in  Figures  33-36  and  are  plotted  as  a  function 
of  azimuth.  For  reference  the  theoretical  effects  of  mutual  coupling  have  also  been  plotted. 

In  general,  the  bearing  errors  follow  the  same  trend  as  the  predicted  mutual  coupling 
error  curves  except  with  a  certain  amount  of  scatter.  There  are,  however,  a  number  of 
exceptions.  In  Figure  34  an  interesting  feature  occurs  between  52°  and  55°  azimuth.  This 
feature  resulted  when  the  transmitter  turned  the  corner  on  Riddell  Road  (see  Figure  4)  and 
began  moving  almost  directly  away  from  the  DF  array  with  a  slightly  decreasing  azimuth 
bearing.  The  result  was  a  repeat  of  some  of  the  azimuth  angles  already  measured  but  for 
different  ranges.  The  resultant  error  in  this  case  appears  to  be  range  dependent.  This  is 
the  kind  of  effect  that  would  be  expected  from  multipath,  particularly  if  the  multipath 
source  produces  very  little  scattering  (i.e.  more  like  a  perfect  reflector  than  an  isotropic 
scatterer). 

In  Figure  35  measurement  108  is  not  shown  due  to  the  large  bearing  error  (—68°). 
During  this  measurement  the  transmitter  was  behind  a  small  building.  This  seems  to  be 
the  clearest  case  of  multipath  given  that  the  direct  signal  path  was  blocked. 

In  Figures  36  there  are  several  places  {<f>  =  —130°,— 43°,  and  54°)  where  the  error 
changes  by  up  to  2°  over  a  narrow  range  of  azimuth  angles.  Again,  the  most  likely  cause 
being  multipath. 

The  amount  of  scatter  observed  for  each  of  the  measurement  sets  can  be  quantified 
by  subtracting  the  predicted  mutual  coupling  errors  from  the  actual  errors.  This  has 
been  done  and  the  RMS  error  results  tabulated  in  Table  7.  Given  the  magnitude  of  the 
RMS  errors,  it  would  appear  that  multipath  y/as  the  dominant  contributor  since  all  other 
known  error  mechanisms  have  essentially  been  ruled  out. 

Finally,  in  comparing  the  November  and  March  error  results,  the  scatter  in  the  March 
errors  is  also  almost  twice  that  of  the  November  errors.  The  reason  for  the  larger  scatter 
could  be  due  to  any  number  of  reasons  which  may  have  effected  the  multipath  conditions. 
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20  25  30  35  40  45  50  55 

Azimuth  Angle  (degrees) 

Figure  33:  Revised  error  results  for  November  5,  1992.  The  dots  represent  actual  error 
and  the  solid  line  represents  predicted  error  due  to  mutual  coupling. 
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15  20  25  30  35  40  45  50 

Azimuth  Angle  (degrees) 

Figure  34:  Revised  error  results  for  November  6,  1992.  The  dots  represent  actual  error 
and  the  solid  line  represents  predicted  error  due  to  mutual  coupling. 
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Figure  35:  Revised  error  results  for  March  24,  1993.  The  dots  represent  actual  error  and 
the  solid  line  represents  predicted  error  due  to  mutual  coupling. 
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-50  0 

Azimuth  Angle  (degrees) 


Figure  36:  Revised  error  results  for  March  25,  1993.  The  dots  represent  actual  error  and 
the  solid  line  represents  predicted  error  due  to  mutual  coupling. 
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This  includes  different  environmental  conditions,  different  frequencies,  different  ranges, 
and  different  azimuth  angles.  In  terms  of  environmental  conditions,  the  ground  state 
may  have  had  an  impact  on  the  results  obtained  during  the  field  trials.  For  example, 
in  November  the  ground  was  wet  while  in  March  it  was  frozen  and  covered  with  snow. 
These  differences  would  be  expected  to  cause  differences  in  the  amplitude  and  phase 
of  the  ground  reflected  signal  and  possible  differences  in  the  measured  bearings  for  the 
same  transmitter  positions.  Changing  the  transmitter  frequencies  would  not  be  expected 
to  change  the  amount  of  scatter  in  the  error  results.  However,  this  possibility  cannot 
be  completely  ruled  out  imtil  further  tests  are  done.  Based  on  the  multipath  models 
developed  in  this  report,  shorter  transmitter  ranges  might  be  expected  to  increaise  the 
scatter  by  increasing  the  relative  contribution  of  multipath  generated  ne3U‘  the  transmitter 
(i.e.  the  loss  function  has  less  effect  at  short  range).  Finally  the  signals  received  during 
the  November  trials  traversed  the  field  where  Osprey  system  was  located  from  the  West. 
In  this  direction,  there  are  relatively  few  prominent  physical  features  for  several  hundred 
meters.  By  comparison,  the  equivalent  distances  for  the  March  trials  were  considerably 
less.  For  example,  at  ^  =  270°  a  ditch  and  row  of  bushes  was  located  within  a  few  meters 
of  the  array. 


6.0  CONCLUSIONS 

This  report  discusses  and  analyzes  the  results  of  field  trials  carried  out  to  investigate 
environmental  effects,  in  particular  multipath,  on  an  experimental  eight-channel  VHF 
DF  system  (Osprey).  The  theoretical  effects  of  error  mechanisms  which  degraded  the  DF 
measurements  were  analyzed  in  detail  including  modelling  of  the  multipath  environment. 
This  analysis  was  then  used  to  breah  down  the  measurement  errors  in  terms  of  these 
mechanisms.  The  RMS  error  results  for  the  March  1993  trials  (which  generally  reflect  the 
worst  case  results)  are  shown  in  Table  8.  These  results  have  been  normalized  to  reflect 
accuracy  for  a  signal  at  0  degrees  azimuth. 

Quzmtitatively,  multipath  propagation  was  found  to  be  the  greatest  source  of  error 
followed  closely  by  mutual  coupling.  Based  on  the  modelling  studies  of  both  phenomena, 
important  qualitative  differences  were  found.  Mutual  coupling  errors  are  a  function  of 
transmitter  azimuth  only  (when  a  plane  wave  signal  is  considered)  and  are  predictable. 
The  error  pattern  changes  relatively  slowly  cycling  once  from  the  maximum  positive  er¬ 
ror  to  the  maximum  negative  error  over  a  range  of  180  degrees  in  azimuth.  Multipath 
errors  are  a  function  of  both  transmitter  azimuth  and  range.  The  error  pattern  becomes 
more  unpredictable  (noise-like)  as  the  range  of  the  multipath  source  from  the  DF  system 
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increases  (at  a  range  of  0  meters  multipath  and  mutual  coupling  can  be  regarded  as  the 
same  phenomena).  These  qualitative  differences  make  it  relatively  easy  to  distinguish 
mutual  coupling  effects  from  multipath  effects. 

In  examining  the  multipath  induced  errors,  two  important  observations  were  made. 
The  first  was  that  the  multipath  errors  were  very  noise-like,  indicating  that  significant 
miiltipath  sources  ranged  up  to  several  hundred  meters  or  more.  The  second  was  there 
appeared  to  be  a  transmitter  range  dependency  even  for  transmitter  ranges  over  5  km. 
Based  on  the  theoretical  models,  this  implies  that  sources  of  multipath  which  were  reflec¬ 
tive  in  nature  (i.e.  redirect  the  incident  signal  in  a  new  direction  with  limited  scattering) 
were  significant. 

These  last  two  observations  lead  to  complications  for  site  calibration  if  improved  DF 
accuracy  is  the  objective.  They  indicate  that  calibration  of  the  array  steering  vector  (or 
array  manifold)  at  a  particular  frequency  needs  to  be  performed  as  a  function  of  both 
range  and  azimuth  ,  not  azimuth  as  is  usually  done.  Further  experimentation  will  be 
required  to  verify  these  observations  before  any  definite  conclusions  can  be  drawn. 

Drift  effects  (i.e.  changes  in  the  gain  and  phase  matching  of  the  DF  receiver  channels) 
were  also  observed  which  were  probably  due  to  the  equipment,  but  may  have  also  been 
due  to  changes  in  the  multipath  environment.  Further  experimentation  will  be  required 
to  verify  whether  the  environment  is  a  factor. 

In  general  the  results  described  in  this  report  indicate  that  the  multipath  environment 
is  very  complex.  Modelling  can  be  used,  with  some  success,  to  better  understand  the  sta¬ 
tistical  effects  of  multipath.  For  future  work,  the  theoretical  models  need  to  be  refined  to 
better  approximate  real  life  effects.  Further  field  testing  will  also  be  required  to  verify  the 
models  and  determine  appropriate  model  parameters  (e.g.  density  of  multipath  sources, 
relative  importance  of  reflective  versus  scattering  types,  etc.).  This  could  ultimately  lead 
to  the  development  of  new  approaches  to  VHF  DF  including  (but  not  restricted  to)  new 
algorithms  and  new  calibration  techniques. 


Table  8:  RMS  Error  Summary 


Mechanism 

RMS  Error  (degrees) 

Transmitter  Position 

0.2 

Noise 

0.1 

Wind 

0.1 

Equipment  Drift 

0.2 

Mutual  Coupling 

0.7 

Multipath 

0.8 
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APPENDIX 


A.O  Stochastic  Calibration  Method 

In  the  stochastic  model,  multipath  is  assumed  to  be  completely  unpredictable  and  can  be 
modelled  as  a  spatially  white  Gaussian  process.  The  array  response  equation  for  this  case 
is  discussed  in  Section  3.2  and  given  by  equation  (10)  which  is  reproduced  here  as, 

A-*X  =  SB  +  N  (A.l) 


Given  that  the  errors  produced  by  multipath  are  represented  by  N,  then  a  least  squares 
solution  for  A  and  B  could  be  determined  by  minimizing  the  total  multipath  error  power 
given  by  trace(N^N).  This  is  not  optimal,  however,  as  the  errors  represented  by  N 
are  not  spatially  white  since  the  power  of  the  multipath  signals  are  proportional  to  the 
signal  power.  Therefore,  a  strong  signal  would  be  expected  to  produce  larger  errors  than 
a  weaker  signal.  To  whiten  the  multipath  noise  it  is  necessary  to  normalize  the  noise 
according  to  the  signal  power,  i.e.,  multiply  the  column  of  N  by  Accordingly, 
equation  (A.l)  becomes 

A-^XB-^  -  S  =  NB-^  =  W  (A.2) 

Therefore  a  more  optimum  estimation  of  A  and  B  can  be  achieved  by  minimizing  the 
error  function, 

MSE  =  trace(W"W)  (A.3) 

A  fundamental  ambiguity  exists  with  the  above  approach.  The  problem  is,  given 
that  A'  and  B'  are  solutions  for  A  and  B  which  minimize  the  MSE,  then  so  are  the 
solutions  «A'  and  k“^B'.  To  remove  this  ambiguity,  the  complex  gain  coeflBcients  are 
chosen  so  that 

imag{ao}  =  0  (A.4) 

and  ^  ^ 

i;  |a.f  =  N  (A,5) 

kzzO 

A  closed  form  solution  for  determining  A  and  B  such  that  the  MSE  is  minimized 
is  unknown.  The  actual  approach  used  to  solve  this  problem  was  to  minimize  the  MSE 
iteratively,  first  solving  for  A  only,  then  solving  B  only,  and  repeating.  The  advantage  of 
this  approach  is  that  the  separate  minimizations  involved  have  closed  form  solutions.  For 
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each  successive  iteration,  ooj  Oi, . . .  ,0^-1  can  be  solved  using, 


]Cn=0  PinPin 
52n=0  ^inPin 


(A.6) 


for  0  <  t  <  iV  —  1  where  aj  is  the  unnormalized  estimate  of  Oj,  p,n  and  s,„  are  elements 
of  the  matrices  P  and  S  respectively,  and 

P  =  XB-^  (A.7) 


The  complex  gain  coefficients  are  then  normalized  so  that 


a.  = 


N 


Similarly  the  parcimeters  60, 61, ... ,  hr-i  can  be  solved  each  iteration  using, 

Y'TV— 1 

L  _  Z^tnsO  "mk^rnk 

Ok  —  , 

^rnk'imk 

for  ifc  =  0, 1, . . . ,  r  —  1  where  qmk  is  an  element  of  the  matrice  Q  which  is  given  by. 


(A.8) 


(A.9) 


Q  =  A-*X 


(A.IO) 


The  minimization  process  is  started  by  initializing  the  diagonal  elements  of  B 
according  to 

hk  =  \/x(<fc)^x(tfc)  (A.ll) 

and  then  for  each  iteration,  updating  the  elements  of  A  followed  by  the  elements  of  B, 
according  to  equations  (A.6)-(A.10).  Ten  iterations  was  found  to  be  more  than  enough 
to  ensure  the  accuracy  of  the  final  estimates  were  sufficient  for  this  report. 
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